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PROGRAM  OBJECTIVE  AND  SUMMARY  OF  RESULTS 


The  intent  was  that  work  performed  under  this  contract  would 
increase  our  understanding  of  the  basic  mechanisms  involved  when 
electrons  interact  with  matter  in  its  condensed  pnases.  Emphasis  has 
been  placed  on  studying  the  interactions  occurring  in  those  materials 
which  are  currently  used,  or  which  are  being  considered  for  use,  in 
electronics  systems.  Knowledge  of  quantities  such  as  the  cross  sections 
for  the  various  types  of  electron  interactions  in  these  materials  and 
the  energy-loss  spectra  for  electrons  generated  by  ionizing  radiations 
incident  on  these  materials,  enable  predictions  to  be  made  concerning 
the  performance  of  electronics  systems  in  various  situations,  such  as 
exposure  to  nuclear  radiatio  .  Our  approach  to  this  general  problem 
has  been  both  theoretical  and  experimental.  In  our  experimental  program, 
we  have  made  observations  which  can  be  related  to  electron  interactions 
in  solids.  The  theoretical  portions  of  the  program  have  included  predic¬ 
tive  calculations  for  a  variety  of  interaction  phenomena  as  well  as 
interpretative  studies  in  support  of  our  experimental  efforts.  Results 
from  the  basic  theoretical  and  experimental  programs  have  been  employed 
in  more  specific,  applied  calculations  with  direct  relevance  to  elec¬ 
tronics  systems'  vulnerability. 

Specifically  we  have:  (a)  made  theoretical  predictions  for  energy- 
loss  spectra  of  electrons  transmitted  through  thin  films  of  carbon  and 
compared  these  predictions  with  available  experimental  data;  (b)  deter¬ 
mined  the  optical  properties  of  polyethylene  in  the  photon  energy  range 
from  0.5  eV  to  76  eV;  (c)  studied,  experimentally  and  theoretically,  the 
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optical  emission  from  several  metals  bombarded  by  low-energy  electrons; 
(d)  produced  a  tabulation  of  inverse  mean  free  paths,  stopping  powers, 
csda  ranges,  and  straggling  in  polystyrene  for  electrons  of  energy  — 

10  keV;  (e)  employed  optical  guided  wave  techniques  to  measure  the 
change  in  refractive  index  of  LiF  as  a  function  of  time  after  deposition 
of  the  film;  (f)  demonstrated  a  new  method  for  directly  observing  and 
recording  the  dispersion  relation  for  optical  guided  waves;  and  (g) 
studied  theoretically  some  extensions  and  modifications  of  classical- 
binary-coll ision  cross  sections.  A  detailed  description  of  these  topics 
in  the  form  of  self-contained  chapters,  makes  up  the  body  of  this  report 


CHAPTER  I 


Straggling  and  Plasmon  Excitation  in  the  Energy  Loss  Spectra 
of  Electrons  Transmitted  through  Carbon* 

J.  C.  Ashley,  V.  E.  Anderson,  J.  J.  Cowan,  and  R.  H.  Ritchie 
Health  and  Safety  Research  Division,  Oak  Ridge  National  Laboratory 

Oak  Ridge,  Tennessee  37830 

and 

J.  Hoelzl 

Department  of  Physics ,  Gesamthochschule  Kassel 
Kassel  35,  Germany 

Abstract 

A  model  is  described  for  the  calculation  of  differential  inverse 
mean  free  paths  (DIMFP's)  for  inelastic  interactions  of  electrons  with 
solids.  The  energy  loss  function,  Im[-l/e(q,u>)],  is  represented  as  a 
sum  of  Drude-type  functions.  The  adjustable  parameters  in  the  energy 
loss  function  are  fixed  by  a  fit  to  the  energy  loss  function  in  the 
optical  limit  (q  -*•  0)  obtained  from  data  on  the  optical  constants  of 
glassy  carbon.  The  extension  of  the  optical  energy  loss  function  to 
arbitrary  values  of  q  in  this  model  leads  to  an  analytical  expression 
for  the  DIMFP.  Energy  loss  spectra  for  low-energy  electrons  (*  2000  eV) 
transmitted  through  carbon  foils  of  various  thicknesses  are  determined 
by  solving  the  transport  equation  using  a  convolution  method  and  the 
DIMFP's  determined  as  described  above.  Comparisons  of  these  calculated 
spectra  with  those  measured  by  Jacobi  show  that  the  model  provides  a 
good  description  of  both  the  broad  straggling  distribution  and  the  fine 
structure  due  to  plasmon  excitation. 

★ 

Research  sponsored  jointly  by  the  Deputy  for  Electronic  Technology, 
Air  Force  Systems  Command,  under  Interagency  Agreement  DOE  No.  40-226-70 
and  the  Division  of  Biomedical  and  Environmental  Research,  U.  S.  Depart¬ 
ment  of  Energy,  under  contract  W-7405-eng-26  with  the  Union  Carbide 
Corporation. 

Computer  Sciences  Division 
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I.  INTRODUCTION 


In  order  to  fully  understand  the  interaction  of  ionizing  radiations 
with  matter,  one  needs  to  be  able  to  measure  or  predict  theoretically  what 
happens  to  both  the  particles  and  the  matter  with  which  they  interact.  As 
a  basic  component  in  solving  this  problem,  we  are  applying  theoretical  models 
to  predict  the  distribution  in  energy  of  electrons  which  are  transmitted 
through  layers  of  material  for  electron  energies  and  layer  thicknesses  for 
which  straggling  is  an  important  factor. 

Characteristic  energy  loss  experiments^  typically  use  electrons  with 
energies  £  20  keV.  At  these  energies,  electron  mean  free  paths  are  >  100  A 
so  electrons  transmitted  through  films  of  a  few  hundred  Angstroms  thickness 
will  undergo  only  a  small  number  of  energy  loss  events.  This  type  of  experi¬ 
ment  allows  one  to  study  the  single  energy  Toss  processes  (both  collective 
and  single-particle)  which  are  characteristic  of  the  given  medium.  For 
lower  electron  energies,  say  £  2000  eV,  mean  free  paths  are  much  smaller 
and  multiple  energy  loss  events  become  important  in  determining  the  spectra 
of  transmitted  electrons  even  in  rather  thin  films.  Such  multiple  energy 

loss  events  lead  to  straggling  as  is  evident  in  the  experimental  work  of 

2 

Jacobi  on  transmission  of  low-energy  electrons  through  thin  carbon  films. 

His  results  for  the  energy  distributions  of  electrons,  Ny(E)  vs  E,  are 
shown  in  Fig.  1  for  monoenergetic,  primary  electron  beams  with  energies 
ranging  from  400  eV  to  1600  eV  transmitted  through  carbon  foils  80  A  and 
230  A  thick.  The  energy  scale  in  the  figure  (horizontal  axis)  may  be 
delineated  by  the  positions  of  the  elastic  scattering  peaks.  The  structure 
at  electron  energies  ^  20  eV  below  the  elastic  peak  results  from  Auger  and 
plasmon  excitations  in  the  material.  At  very  low  energies,  secondary 
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electrons  also  contribute  to  the  distributions.  The  above  contributions 

are  also  seen  in  the  energy  distributions  of  reflected  electrons.  The 

broad  asymmetric  component,  which  appears  only  in  transmission,  has  the 

characteristic  of  becoming  narrower  and  shifting  to  higher  energies  as 

either  the  foil  thickness  is  decreased  or  the  Drimary  energy  is  increased. 

This  component  is  oue  to  straggling  of  the  primary  electrons  as  they  are 

degraded  in  energy  in  passing  through  the  foil. 

Employing  the  simple  Bohr  impact  parameter  model  for  electron- 

electron  interactions  and  the  straggling  theory  of  Landau,  the  strag- 

3 

gling  distributions  could  be  accounted  for  in  a  reasonable  manner. 

However,  the  structure  due  to  multiple  plasmon  excitations  superimposed 
on  the  broad  straggling  distributions  was  not  accounted  for  in  the 
simple  model  used  in  reference  3.  Since  the  experimental  energy  spectra 
indicate  substantial  contributions  from  plasmon  excitation,  a  more 
realistic  description  of  the  interaction  of  the  incident  electrons  with 
the  foil  should  include  the  plasmon  excitation  process  as  well  as  single 
particle  collisions. 

In  the  work  presented  here  we  will  concentrate  on  the  electron 
energy  loss  region  where  plasmon  excitations  are  important  in  deter¬ 
mining  the  structure  of  the  energy  1oss  spectrum  (electron  energy  losses 
up  to  ~ 200  eV).  The  theory  and  models  used  in  the  calculations  will  be 

described  in  the  next  section.  The  theoretical  predictions  will  be 

2 

compared  with  the  experimental  results  for  carbon  in  Section  III 
followed  by  some  general  observations  on  the  validity  and  usefulness  of 
these  calculations  for  predictions  in  applied  problems. 
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II.  THEORETICAL  DESCRIPTION  OF  ELECTRON  ENERGY  LOSSES 


A  basic  component  of  the  theory  is  a  suitable  description  of  the 
differential  inverse  mean  free  path  (DIMFP),  d'  / dnj ,  for  inelastic 
interactions  of  an  electron  with  the  solid.  Since  we  are  interested  in 
the  distribution  of  energy  losses  only,  the  DIMFP,  or  the  .  oDqo  litv 

4 

per  unit  path  length  per  unit  energy  for  an  energy  loss  j,  is  given  by 


f?1 


cLA'/dcu  h  t(E,co):=  j  ^  ImC'i/sC^cj) 0) 

9 

where  E  "  v  /2  is  the  energy  of  the  electron,  v  is  the  electron  velocity, 

and  the  1’mits  on  the  integration  over  momentum  transfer  q  are  given  by 

— 

q =  v  •  \v  -  2 if.  In  the  discussion  here  we  work  in  atomic  units 
(a.u.)  where  fi  =  m  =  e  =  1;  conversions  to  other  units  will  be  given 
where  appropriate.  The  energy  loss  function  Im[ - 1 /e (q,u ) ] ,  the  imaginary 
part  of  the  negative  reciprocal  of  the  complex  dielectric  response 
function,  c,  of  the  solid,  describes  the  response  of  the  medium  to 
energy  and  momentum  transfers.  A  simple  model  for  the  dielectric 
response  function  which  emphasizes  the  influence  of  single-particle 
interac’ions  for  large  momentum  transfer  and  interactions  with  plasmons 
for  small  momentum  transfer  is  given  by  a  sum  of  Drude-type  functions  as 


li n  =  2 


A  LV[cc) 


(2) 


An  energy  loss  function  of  this  form  was  used  by  Hamm  et  al.  for  calcu¬ 
lations  of  electron  DIMFP's  in  liquid  water.  In  the  limit  q  -*•  0  (optical 
limit),  the  constants  A^ ,  ,  and  may  be  obtained  from  optical  data 

for  the  medium  or  from  characteristic  energy  loss  dataJ’^ 
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Using  Eq.  (2),  where  q  f  0,  the  integration  over  momentum  transfer 
to  obtain  the  DIMFP  may  be  performed  analytically  to  give 


(3) 


where 

7"  (cL)  =  L  a ScOol «<.  +  -  woc  fr  ( (J0ZL  -  cj2)]  L  L 

+  K(<JoL*L+  CtX U)Z  - O)2)  -t-  CLlc0^l  ]  B; 

with  the  definitions 


<2-6  =  ^  ; 

C;  =  (cvH-t-ozylH , 

ft=LiO-«Vcnife, 

xi*.  C(x?  +yf)  /  Cxf.+^o  , 

B^WVy/x^)  -  ta.rCl(fr/xiJ  , 

yi  =  -c;A  ,  +  -ci«i  • 

In  the  equation  for  ,  the  angles  e|  h  tan"  (y^/x^)  must  lie  in 
.  + 

the  region  -it  <  e .  <_  tt. 

Having  obtained  values  for  t(E,w),  one  can  calculate  the  energy 
loss  distribution  f(t,w)  which  gives  the  probability  that  an  electron  of 
energy  E  will  have  lost  energy  w  on  traversing  a  thickness  t  of  the 
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solid.  The  energy  loss  distribution  is  a  solution  of  the  transport 
equation 


df  (x;dj)/dx  -J  v(EjCo')  fUjid-co^ctd)'  -ffyuj) A'(E)  ) 


where  A’  (E)  is  the  inverse  mean  free  path  given  by 


A'CE)  =  I  T(E,oo)du>  .  (5) 

Jo 

Values  of  f(t,w)  were  obtained  using  a  convolution  technique  developed 
by  Williams.7  This  technique  has  been  employed  recently  in  a  detailed 
study  of  straggling  distributions  for  protons  in  aluminum.  The  solution 

g 

of  the  transport  equation  may  be  written  in  the  form 


that  is,  the  distribution  after  traversing  a  thickness  2x  of  the  solid 
is  obtained  as  a  convolution  of  the  distribution  of  electrons  which  have 
lost  energy  w  -  w'  in  thickness  x  with  the  distribution  for  energy  loss 
ui*  in  x.  We  start  with  a  thickness  xQ  for  which  the  probability  of  a 
single  energy  loss  event  is  quite  small.  The  distribution  for  this 
situation  is  taken  to  be 


0-  x0A')S(co)  +  x0rfE,oj)  . 


(7) 


In  this  distribution  the  zero  energy  loss  term  [given  by  the  Dirac  delta 
function  <5(w)]  has  been  decreased  by  the  probability  for  any  energy  loss 
event,  x  A  \  That  energy  loss  event  was  assigned  to  the  energy  loss 
distribution  as  determined  by  the  DIMFP  through  the  term  xqt.  Given 
this  initial  distribution,  Eq.  (7),  f(t,w)  may  be  obtained  for  any 
desired  value  of  t  by  repeated  application  of  Eq.  (6).  The  result  after 
n  convolutions  may  be  written 

f(t,C j)=  bZ  8( 0))  4-  Fn(oj)  ;  (8) 

where 

b=  l-  xqA'  , 

Fn(w)=2 bz  Fn_,(uj)  +f  F^Coj')  du)' 

o 

and 

XoT(£;Cj)  . 

The  energy  loss  distributions  calculated  using  this  procedure  are  expected 
to  be  valid  for  energy  losses  small  compared  to  the  incident  electron 
energy  E  since  the  DIMFP  determined  for  E  is  used  in  the  convolutions. 

III.  CALCULATIONS  OF  ENERGY  LOSS  SPECTRA 

To  evaluate  the  constants  required  in  Eq.  (2),  we  used  the  optical 

g 

data  of  Williams  and  Arakawa  on  glassy  carbon  for  photon  energies  up  to 
80  eV.  The  energy  loss  function  determined  from  this  data  is  shown  as 
the  solid  curve  in  Fig.  2  and  corresponds  to  carbon  of  density  p  =  1.47 
g/cm  .  The  energy  loss  in  electron  volts  is  given  by  27.2  x  w  (a.u.). 
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<  • 
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The  peaks  in  this  function  are  well  understood  in  terms  of  collective 
oscillations  of  the  a  and  v  electrons.  Using  Eq.  (2)  for  q  =  0,  we 
determined  values  of  the  parameters  A- ,  7^,  and  u>  .  using  a  sum  of  four 
functions  and  fitting  to  the  optical  energy  loss  function.  The  result 
is  shown  as  the  dashed  curve  on  Fig.  2.  At  energies  below  the  K-edge  in 
carbon  at  ~  282  eV,  the  optical  energy  loss  function  is  determined  by 
the  response  of  the  valence  electrons.  As  a  check  on  the  fitting  pro¬ 
cedure,  we  require  that  Eq.  (2)  account  for  all  the  valence  electrons 
through  the  sum  rule 

f  cj  ImL-i /stored)] duj  =■  Z7 rznv  J  (9) 

where  ny  is  the  valence  electron  ^density  assuming  a  contribution  of  four 
electrons  per  carbon  atom.  From  Eq.  (9)  with  Eq.  (2),  we  have  £  A^  = 
4nnv.  Energy  losses  due  to  K-shell  ionization  are  not  included  since 
they  are  outside  the  energy  loss  range  studied  here. 

Solid  carbon  occurs  in  several  forms  with  a  wide  range  of  densities. 

2 

The  evaporated  films  used  in  the  work  of  Jacobi  correspond  to  a  density 
3 

P  =  2.0  g/cm  .  Ideally,  we  would  use  optical  data  obtained  on  the  same 
form  and  density  of  carbon  used  in  this  experiment  to  determine  the 
parameters  in  Eq.  (2).  Since  such  information  is  not  available,  the 
energy  loss  function,  Eq.  (2),  determined  from  glassy  carbon  data  must 
be  modified  to  correspond  to  this  higher  density.  From  characteristic 
energy  loss  experiments,^  it  is  found  that  the  peaks  due  to  excitation 
of  collective  oscillations  or  plasmons,  as  seen  in  Fig.  2  for  qlassy 
carbon,  are  predicted  quite  well  for  different  densities  assuming  the 
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position  is  proportional  to  /p.  The  values  of  <j  ■  have  each  been  scaled 

0 1  /I 

to  p  =  2  g/cnr  using  this  proportionality.  Since,  from  Eq.  (9),  ^ 

i  =  l 

4irnv<  we  have  scaled  each  A.  by  the  factor  2/1.47.  The  values  of 

which  are  related  to  the  widths  of  the  peaks,  are  assumed  to  be  proportional 

to  ,  as  indicated  by  arguments  based  on  a  free-electron  approximation.^ 

The  set  of  parameters  determined  by  the  fit  to  the  glassy  carbon  data. 


Fig.  2,  and 

scaled  to  a  density  , 

=  2.0  g/cnr 

are  shown  in  Table  I, 

expressed  in 

atomic  units. 

Table  I. 

Parameters  Used  in 

Calculations 

of  Energy  Loss  Spectra 

i 

Ai_ _ 

Ui  . 

oi _ 

_ ^i_ _ 

1 

0.00736 

0.245 

0.100 

2 

0.3755 

0.9005 

0.400 

3 

0.2796 

1 .458 

0.9007 

4 

0.08463 

2.573 

1.00 

Examples  of  DIMFP's  calculated  from  Eq.  (3)  using  the  parameters 
given  in  Table  I  are  shown  in  Fig.  3,  with  x  and  w  in  atomic  units,  for 
electrons  of  energy  500,  1000,  and  2000  eV.  The  structure  in  the  energy 
loss  function,  Fig.  2,  is  also  present  in  the  D IMF P  with  maxima  at~0.25 
(6.8  eV)  and  at  ~0.93  (25.3  eV)  corresponding  to  the  plasmon  losses  in 

3 

carbon  at  the  density  p  =  2.0  g/cm  . 

Electron  inelastic  mean  free  paths  have  been  calculated  from  Eqs. 

(5)  and  (3)  using  the  parameters  in  Table  I.  The  results  are  given  by 
the  solid  line  in  Fig.  4  where  we  plot  A  in  Angstroms  as  a  function  of 
electron  energy  in  electron  volts.  The  dashed  curve  was  calculated  from 
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1  ?  3 

Penn's  electron  gas  model  "  for  carbon  with  c  =  ?.  g/cm  ;  Penn's  values 
fall  about  AO"!  below  ours.  The  experimental  values  of  electron  atten¬ 
uation  lengths  shown  in  the  figure  are  7.5  t  1.5  A  at  263  eV  from  Jacobi 
1  3 

ard  Hoelzl  and  the  points  at  920  eV  and  1169  eV  are  from  Eteinha  dt, 

1 4 

Hudis,  and  Perlman;  no  error  limits  were  assigned  to  the  latter  va1 
The  DIMFP's  predicted  by  this  model  were  used  to  calculate  the 
distribution  of  energy  losses  for  electrons  transmitted  tr, rough  carbon 
foils  of  various  thicknesses  using  the  convolutic  method  described  • 

Section  II.  In  npp  following  figures  we  show  the  ener j  lo  >  di  tr -'but  i-on 
t ( t ,m)  as  3  function  of  energy  loss  for  electrons  o*  energy  E  trans¬ 
mitted  through  a  carton  foil  of  thickness  t.  In  Figs.  5  and  6  w a  can 
see  the  changes  in  the  energy  loss  distribution  for  L  =  1000  eV  as  the 

thickness  of  the  carbon  layer  increases  from  10  A  through  115  A.  The 

2n  f 

contribution  to  f(t,«)  from  zero  energy-loss  electrons,  b  f  5(c)du,  is 
given  approximately  by  e  Even  for  the  small  thickness  t  =  10  A, 

over  30  of  the  electrons  have  lost  energy  on  transmission  through  the 
carbon.  In  cig.  5,  the  plasmor,  loss  peak  at -0.93  at  first  grows  with 
increasing  t  but  then  decreases  as  maxima  begin  to  develop  at  multiples 
of  this  loss.  For  a  thickness  of  20  A,  slightly  less  than  one  electron 
mean  free  path  (A  -  26  A  at  E  -  1000  eV),  a  shoulder  corresponding  to 
electrons  that  have  undergone  two  losses  is  seen  at  w  =?  1.9.  Multiples 
of  the  0.25  energy  loss  are  too  small  to  be  seen  in  these  curves.  Mul¬ 
tiple  plasmon  losses  become  more  evident  for  the  larger  values  of  t 
shown  in  Fig.  6.  The  growth  of  the  straggling  distribution,  with  struc¬ 
ture  associated  with  multiple  plasmor,  excitations,  is  clearly  shown  in 
these  figures. 
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It  is  important  to  point  out  that  our  theoretical  model  assumes  the 
thin  films  are  uniform  and  homogeneous,  and  the  dielectric  properties  of 
the  carbon  are  the  same  as  for  the  bulk  material  of  a  given  density. 
Figures  5  and  6  illustrate  the  development  of  the  energy  loss  distribu¬ 
tion  based  on  these  assumptions.  Experimentally,  it  is  quite  difficult 
to  obtain  a  uniform  film  with  a  thickness  a  100  A;  this  problem  should 
be  kept  in  mind  wher1  making  comparisons  of  theory  with  experiment. 

In  Fig.  7  we  have  compared  our  calculations  with  some  representative 
experimental  results  of  Jacobi  (dashed  curves)  with  80-fl  carbon  foils 
for  incident  electron  energies  from  600  eV  tnrough  1600  eV.  Data  taken 
from  Fiq.  1  on  the  number  of  transmitted  electrons  of  1600-eV  incident 
energy  were  scaled  by  a  factor  which  made  the  overall  maximum  roughly 
the  same  height  as  predicted  by  the  theory.  The  same  factor  was  used  to 
scale  the  transmitted  intensities  for  the  lower  energies.  Electron 
energy  loss  was  measured  from  the  elastic  peak.  As  one  goes  toward 
zero  energy  loss  in  Fig.  7,  the  experimental  data  begin  to  rise  into 
a  slightly-broadened  elastic  peak  while  the  theoretical  calculations 
approach  zero  [since  the  theoretical  model  assumes  the  elastically 
scattered  electrons  appear  at  exactly  zero  energy  loss  as  given  by 
the  delta-function  term  in  Eq.  (8)].  Remarkably  good  agree¬ 
ment  is  found  between  theory  and  experiment  for  the  1600- ,  1400-,  1 200- , 
and  1000-eV  comparisons.  The  positions  of  the  plasmon  peaks  and  the 
overall  shapes  of  the  data  are  predicted  quite  well  by  the  theory.  As 
one  goes  to  lower  incident  electron  energies,  the  positions  of  the 
maxima  in  the  broad  straggling  distributions  are  predicted  quite  closely 
by  the  theory  down  through  600  eV.  At  400  eV  (not  shown  in  Fig.  7),  the 
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i.  os  i  t  on  of  the  maximum  is  predicted  to  occur  at  — 1 93  eV ,  while  the 
experimental  data  give  a  value  of  — 1 60  eV.  This  represents  quite  an 
improvement  over  the  results  of  reference  3  where  the  peak  positions 
predicted  from  a  simple  theoretical  model  were  about  twice  as  large  as 
those  iound  experimentally  for  the  80-ft  foil.  Our  calculations  for  a 
2 30- A  Ccroon  foil  predict  the  position  of  the  maximum  in  the  energy  loss 
spectrum,  or  the  most  probab1e  energy  loss,  to  occur  at  a  va^je  ~1.8 
times  that  found  experimentally  (see  Fig.  1)  for  incident  electron 
energies  of  1000,  1200,  1400,  and  1600  eV. 

The  calculated  f(t,u>)  tend  to  overestimate  the  experimental  intensities 
at  the  lower  electron  energies.  The  ra^io  of  theoretical  to  experimental 
values  (at  the  maxima)  is  1.3  at  800  eV,  1.6  at  600  eV,  and  3.8  at  400  eV 
for  the  80-A  foil  (see  Fig.  7).  As  pointed  out  earlier,  the  theory  used 
for  calculation  of  the  straggling  distributions  assumes  that  energy  losses 
are  small  compared  to  the  incident  electron  energy.  This  condition  is 
most  nearly  fulfilled  for  the  highest  energy  (1600  eV)  for  the  range  of 
energy  losses  shown  in  Fig.  7  (  from  zero  up  to  —  163  eV).  For  example, 
a  1600-eV  electron  which  has  lost  100  eV  has  its  mean  free  path  reduced 
by  6  while  this  energy  loss  reduces  the  mean  free  path  of  a  600-eV 
electron  by  13, .  This  results  in  a  larger  increase  in  the  probability 
for  a  further  energy  loss  event  to  occur  for  the  electron  of  initial 
energy  600  eV  than  for  the  electron  of  initial  energy  1600  eV.  This 
would  account  for  the  increase  in  the  ratio  of  the  theoretically  pre¬ 
dicted  maximum  to  that  observed  experimentally  as  one  goes  to  lower 
incident  electron  energies,  since  the  theory  uses  a  constant  value  of 
electron  mean  free  path  corresponding  to  the  incident  electron  energy. 


The  analytic  form  of  the  DIMFP  described  in  this  work  should  be 
useful  in  conjunction  with  a  Monte  Carlo  solution  of  the  transport 
equation  to  account  for  the  change  in  DIMFP  after  each  energy  loss  event 
experienced  by  an  electron.  Using  this  approach  we  could  also  study  the 
build-up  of  secondary  and  Auger  electrons  and  their  contributions  to 
the  energy  loss  spectra.  We  plan  to  investigate  this  approach  for 
calculation  of  energy  loss  distributions  for  low-energy  electrons. 

IV.  SUMMARY 

A  model  has  been  described  for  the  calculation  of  differential 
inverse  mean  free  paths  of  electrons  interacting  with  a  medium.  Optical 
data  for  glassy  carbon  was  used  to  determine  the  parameters  in  the  model. 
The  DIMFP's  thus  determined  were  employed  in  calculations  of  energy  loss 
distributions  for  electrons  transmitted  through  thin  layers  of  carbon. 
Comparisons  with  available  experimental  data  indicated  that  the  theoretical 
description  of  the  energy  loss  spectra  accounts  quite  well  for  the 
structure  in  the  straggling  distributions  due  to  plasmon  excitation. 
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FIGURE  CAP' IONS 


Fig.  1.  Energy  distributions  of  electrons  transmitted  through  carbon 

foils  (thicknesses  80  A  and  230  A)  for  primary  electron  energies 
from  400  to  1600  eV.  These  curves  are  from  the  work  of  Jacobi, 
reference  2. 

Fig.  2.  The  optical  energy  loss  function  for  glassy  carbon  (solid  curve) 
as  derived  from  experimental  data,  reference  9,  vs  photon  energy 
in  atomic  units  (bottom  scale)  or  electron  volts  (top  scale). 

The  dashed  curve  shows  our  fit  to  the  experimental  result  using 
a  sum  of  Drude-type  terms. 

Fig.  3.  Differential  inverse  mean  free  paths  as  a  function  of  energy 
loss  for  electrons  of  energy  500,  1000,  and  2000  eV.  These 
curves  are  based  on  the  parameters  obtained  from  the  theoretical 
fit  shown  in  Fig.  2  (see  text  for  details). 

Fig.  4.  Electron  mean  f_ee  paths,  A  in  Angstroms,  as  a  function  of 
electron  energy  for  carbon.  The  solid  curve  is  the  result 
of  this  work,  the  dashed  curve  is  based  on  an  electron  gas 
model  (reference  12),  and  experimental  va^es  of  electron 
attenuation  lengths  are  from  references  13  !•)  and  14  (c). 

Fig.  5.  Calculated  energy  loss  spectra  for  1000-eV  electrons  transmitted 
through  carbon  layers  of  thickness  10,  20,  and  40  A. 

Fig.  6.  Calculated  energy  loss  spectra  for  1000-eV  electrons  transmitted 
through  carbon  layers  of  thickness  57.5,  80,  and  115  A. 

Fig.  7.  Comparisons  of  calculated  electron  energy  loss  spectra  (solid 
curves)  with  experimental  data  from  reference  2  (dashed  curves) 
for  600-,  800-,  1000-,  1 200- ,  1 400- ,  and  1600-eV  electrons 
transmitted  through  80  A  of  carbon. 
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ABSTRACT 


Emission  spectra  from  Ag  films  irradiated  by  low  energy  electrons 
(20-1500  eV)  have  been  measured,  and  the  results  compared  with  theory. 
For  relatively  smooth  films,  two  peaks  in  the  spectra  are  resolved. 

One  at  3.73  eV,  the  volume  plasmon  energy,  is  attributed  to  transition 
radiation  and/or  bremsstrahlung .  The  second,  at  about  3.60  eV,  is 
very  sensitive  to  surface  roughness  in  both  position  and  magnitude 
and  is  produced  by  roughness -coupled  radiation  from  surface  plasmons. 
For  rough  fi lms ,  the  roughness -coupled  radiation  dominates  the  emission. 
In  addition  .o  spectral  shapes,  the  polarization  of  the  radiation 
and  its  intensity  as  a  function  of  electron  energy  were  measured. 

The  experimental  results  are  compared  with  new  calculations  of 
roughness -coupled  emission  which  account  for  most  of  our  observa¬ 
tions.  They  indicate  that  high  wavevector  roughness  components 
play  the  dominant  role  in  the  emission  process. 


Introduction 

Energetic  electrons  incident  on  a  metal  surface  produce  radiation 
via  several  physical  processes.  Transition  radiation  is  produced  by 
the  collapsing  dipole  field  generated  as  the  particle  enters  the 
surface,  and  bremsstrahlung  is  produced  by  electron  scattering 

in  the  metal. (3‘5)  These  processes  produce  spectra  characterized  by 
a  peak  or  shoulder  at  the  volume  plasma  frequency  (my)  of  the  electron 
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gas  with  a  continuum  of  radiation  extending  to  lower  frequencies. 
They  may  be  distinguished  experimentally  since  transition  radiation 
produces  P-polarized  light  only  (electric  vector  in  plane  defined 
by  the  incident  beam  and  the  surface  normal),  while  bremsstrahlung 
produces  light  with  both  P  and  S  polarizations.  The  electrons  also 
excite  surface  plasma  waves  that  do  not  radiate  so  long  as  the  sur¬ 
face  is  smooth. A  rough  surface  allows  the  plasmons  to  de-excite 
roughness- coupled  surface'plasmon  radiation  (RCSPR) . 1 

Other  experiments  provide  information  about 
one  or  the  other  of  the  two  physical  processes  giving  rise  to  RCSPR 
from  Ag.  Energy  losses  suffered  by  electrons  scattered  from  the 
metal  surface  provide  information  about  the  plasmon  production 
process, while  ref lectivity , C1 3" 1 7)  light  scattering^ 8~20)  and 
absorption ^  experiments  yield  information  about  the  roughness 
coupling  of  plasmons  to  photons. 

Radiation  produced  by  high  energy  electrons  (>  20  keV)  incident 
on  Ag  has  been  studied  by  several  workers . (22~27)  in  this  energy 
range,  transition  radiation  dominates  the  spectra  observed  by  most 
workers , (22-26)  although  RCSPR  has  been  observed  for  sufficiently 
rough  films. (25~27)  in  this  paper  we  report  measurements  of  the 
radiation  spectra  produced  by  low-energy  electrons  (20-1500  eV) 
incident  on  Ag.  For  sufficiently  smooth  surfaces,  both  volume- 
plasmon  radiation  (transition  radiation  and/or  bremsstrahlung)  and 
RCSPR  were  observed,  but  for  most  surfaces  the  RCSPR  dominates  the 
observed  spectra. 
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The  formal  aspects  of  tne  theory  we  use  to  describe  our  results 
are  given  in  a  separate  paper  that  precedes  this  or.e.i?t:'  We  will 
describe  here  some  of  the  physical  processes  contributing  to  the 
formal  results. 

A  dispersion  curve  u  (k)  for  surface  plasinons  is  shown  in 

Fig.  !a  along  with  the  dispersion  line  of  light.  The  plujmon 

curve  approaches  the  light  line  in  the  retarding  region  (k'0.5) 

where  the  plasrr.on  pnase  velocities  become  comparable  with  the 

elocitv  of  light  in  vacuum,  and  approaches  a  limiting  value  w 

in  the  normal  plasmon  region  (ki2.  .  Wavevectors  k  in  the  figure 

ur'  given  in  units  of  —  =  t —  (=  l.S  x  10  '  A  1  for  Ag .  .  For 
c  A 

s 

a  free-electron  gas  with  no  damping,  RCSPR  would  have  a  sharp  upper  limit 
at  An  important  feature  of  surface  plasmons  for  a  real  metal  however. 

s  • 

is  that  they  are  strongly  lifetime  broadened,  particularly  at  high  h-vectors  where 
u  (k)  apprcac  :ter  os. C29)  We  will  cite  the  presence  of  large  amounts 
of  RCSPR  with  <j'»u>s  as  t  idence  of  the  importance  of  high  k-vector 
plasmons  to  the  observed  spectra. 

The  surface  roughness  is  described  in  terms  of  an  autocorrelation 
function  which  may  be  fourier  analyzed  into  components  of  wavevector 
Each  q  represents  a  periodic  structure  that  can  couple 
plasmons  and  photons  of  the  same  frequency  when 

£?.l  1  ♦  4  =  £pl»  (1) 

where  kgn  *s  t*ie  component,  of  the  wavevector  of  light  lying  in  the 
radiating  surface.  The  coupling  is  shown  schematically  in  Figure  lb. 

For  an  observation  angle  8  from  the  normal,  each  roughness  component 
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q  can  couple  plasmons  with  wavevector  magnitudes  a-k,  sin  9<k  <q+k 

D 1  X" 

sin  9  to  the  radiation  field.  For  high  a-vector  roughness  components 
(q>>k,  )  or  photon  emitted  near  normal  (sin  0*10 ) ,  kpj^q  so  that  the 
oh served  RCSPR  spectrum  is  determined  in  large  measure  by  the  distri¬ 
bution  and  magnitude  of  the  surface  roughness  components. 

The  other  important  factor  determining  the  RCSPR  spectrum  is  the 
distribution  of  plasmons  in  2-D  k-space  produced  by  an  electron  of 
energy  E  incident  at  an  angle  f  to  the  surface  normal.  The 
formalism  of  Kretschmann  et  al  in  the  accompanying  paper  is 
valid  for  all  regions  of  the  dispersion  curves  but  a  more  trans¬ 
parent  result  an  be  obtained  from  an  earlier  treatment  of  Wilems 
and  Ritchie  which  is  valid  for  the  normal  plasmon  region  where 
mfk)'Vus  at  all  k-vectors.  ^  For  normally  incident  electrons  of 
energy  E,  the  probability  of  producing  a  plasmon  of  wavevector  k  may 
be  written  in  the  form 


P(k)  =  const 


Vq/c ‘ 
(Vo'/cT 


k2 

►  1/k2) 2 


= const! ' 


(2E/mc2 


where  v0  is  the  electron  velocity,  m  the  electron  mass  and  k  is 
written  in  units  of  — .  For  a  given  electron  energy  P(k)  is  pro¬ 
portional  to  k  for  k<< (mc2/2E) 1 '  2 ,  reaches  a  peak  at  k  =  (mc2/2E)1/2 
and  decreases  for  higher  values  of  k.  Since  (mc2/2E)'/2  =  16  for 
E  =  1000  eV,  relatively  few  small  k  plasmons  are  produced  by  electrons 
with  the  energies  used  in  this  study. 

We  find  that  our  experimental  data  can  be  accounted  for  quite 
satisfactorily  if  we  assume  that  the  RCSPR  spectra  are  produced 

primarily  by  a  relatively  narrow  distribution  of  high  q-components  in  the  roughness 


spectraii.  The  well  known  presence  of  another  peak  centered  at  q  =  0 
has  little  effect  on  our  results  since  few  small  k  plasmons  are 
generated  by  low  energy  electrons.  Assuming  a  high  q  roughness  peak, 
equation  (1)  shows  that  the  peak  of  the  RCSPR  spectrum  observed  near 
normal  results  froi:  the  decay  of  plasmons  with  wave- 
vector  kj. .  j  y  q(i  ,  wher  e  q^  is  the  dominant  roughness  component  Thus 
the  peal  intensity  of  the  RCSPR  spectrun  is  proportional 
to  the  probability  of  producing  a  plasmon  of  wavevector 
k  =  o,  .  By  fitting  equation  (2) 

to  a  plot  of  peak  intensity  versus  electron  energy,  it  is  possible 
to  obtain  an  estimate  of  q^,  the  dominant  component  of  the  high  q 
roughness  peak.  The  theory  of  Kretschinann  et_  al .  in  the  accompanying 
paper  develones  the  above  ideas  in  a  consistent  formalism  which 
allows  us  to  obtain  an  estimate  of  the  dominant  high  q  roughness 
component  of  a  Ag  film  from  the  experimental  data. 

\n  important  extension  of  the  present  theory,  however,  is  the 
inclusion  of  the  effect  of  roughness  in  coupling  plasmons  to  each 
other  as  well  as  to  the  photon  field. 3>30)  For  a  single  roughness 
component,  the  calculations  predict  a  prominent  splitting  of  the 
RSCPR  peak  at  The  splitting  is  not  resolved  in  the 

RCSPR  spectra  of  Ag  reported  here,  presumably  due  to  the  finite 
width  of  the  roughness  peak  and  to  the  lifetime  broadening  of  the 
plasmons,  but  it  does  account  in  a  natural  way  for  a  shift  in  the 
position  of  the  peak  of  the  RCSPR  spectrum  that  is  observed  as  the 

(■3 

surface  is  roughened.  A  double -peaked  RCSPR  spectrum  of  potassium  1 
has  been  observed  in  our  laboratory  that  is  consistent  with  the 
prediction  of  the  Kretschmann  theory. 


Experimental  Techniques  and  Results 

A  schematic  diagram  of  the  apparatus  used  in  the  studies  is 
shown  in  Fig.  2.  The  angle  between  the  electron  beam  and  mono¬ 
chromator  port  is  fixed  at  35°.  Most  of  the  results  reported  here 
were  obtained  with  near-normal  incidence  electrons  and  an  observation 
angle  of <40°.  Silver  films  were  prepared  in  situ  by  evaporating 
Ag  onto  various  substrates  with  the  film  thickness  monitored  by  a 
quartz  crystal  thickness  monitor.  Photons  emitted  by  the  films 
were  analyzed  in  the  wavelength  range  between  3000  A  and  5500  A 
using  a  grating  monochromator  and  photomultiplier.  The  exper¬ 
iments  were  performed  in  an  unbaked,  ion-pumped  vacuum  system 
operated  with  a  base  pressure  of  about  5  x  10  7  torr  which  rose 
to  the  10'6  torr  range  during  evaporation. 

Film  preparation  and  treatment  had  a  very  strong  effect  on 
the  radiation  produced  by  the  electrons.  Substantially  different 
results  were  obtained  for  Ag  films  deposited  on  shiny  Ag  sheet, 
shiny  Ta  plate,  Mo  plate  with  a  matte  finish,  and  glass  slides. 

The  rate  of  Ag  deposition  and  local  heating  of  the  electron  beam 
also  affected  the  radiation  observed. 
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Fig.  3  shows  the  emission  spectra  from  a  specularly  reflecting  Ag 


plate  irradiated  with  electrons  of  500  eV  and  1000  eV  both  before 
(dashed  curve)  and  after  (solid  curve)  the  rapid  (M000  A/min) 
evaporation  of  a  2000  A-thick  Ag  layer,  a  double-peaked  spectrum 
is  observed  with  peaks  located  at  approximately  3320  A  (3.73  eV) 
and  3440  A  (3.60  eV)  .  N'ote  that  the  magnitude  of  the  3440  A  peak 
is  enhanced  for  the  evaporated  film.  The  positions  of  the  volume  plasmon 

O 

oscillation  (fiwv  =  3.76  eV,  \y  =  3300  A)  and  of  the  limiting  surface 

O 

plasmon  excitation  (hcus  =  3.66  eV,  As  =  3390  A)  are  indicated  by  arrows 
on  this  and  most  of  the  following  figures. 

Fig.  4  shows  results  from  a  Ag  film  evaporated  on  a  matte 
Mo  plate.  The  double-peaked  spectrum  is  retained,  with  one  peak 

O 

located  at  about  3300  A  as  before,  but  the  dominant  peak  is  now 

o 

centered  at  about  3500  A. 

Fig .  5  shows  the  results  obtained  with  300  eV  electrons 
on  a  rapidly  evaporated  thick  film  of  Ag  on  Ag  sheet  along  with 
spectra  obtained  when  additional  Ag  is  evaporated  very  slowly 

O 

(^100  A/min)  onto  the  surface.  The  results  on  the  original  film 
can  best  be  interpreted  as  two  poorly  resolved  peaks  of  about 

O  O 

equal  magnitude  located  at  3320  A  and  3440  A.  During  the  slow 

O 

evaporation,  the  intensity  at  3320  A  remains  virtually  unchanged,  but  the 

O 

3440  A  peak  grows  dramatically  in  magnitude  and  shifts  to  slightly 
longer  wavelengths. 

Fig  .  6  shows  the  emission  from  Ag  on  glass  for  films  of 

O 

increasing  thickness.  A  rather  broad  peak  centered  at  3440  A 


grows  in  magnitude  and  shifts  to  longer  wavelength  with  increasing 
film  thickness.  In  these  curves  the  sur face-p lasmon  radiation 
completely  dominates  the  spectra,  and  the  short-wavelength  peak  is 
barely  visible  as  a  shoulder. 

If  excessive  electron  beam  power  was  dissipated  in  the  \g  on 
glass  films,  the  appearance  of  the  films  changed  with  time  from 
shiny  to  frosty  as  a  result  of  film  damage  by  local  heating.  The 
frosty  spots  on  such  films  were  rough  whet,  viewed  through  an  optical 
microscope.  Fig  .  7  shows  the  accompanying  change  in  the  emission 
as  a  Ag  film  was  subjected  to  a  high-intensity  (M  mA)  electron 
beam.  The  peak  is  shifted  to  longer  wavelengths  as  film  roughening 
proceeds.  A  fresh  evaporation  of  Ag  on  such  i  damaged  region  does 
not  move  the  emission  peak  back  to  its  original  posit i  n  as  would 

-n 

be  the  case  if  the  shift  were  due  to  the  formation  of  an  cxide  or 
other  contaminant  layer  on  the  surface. 

In  nearly  all  curves  of  Figs.  3  through  6,  there  is  evidence 
for  a  high-energy  component  at  about 

3300  X  which  changes  little  with  changes  in  the  surface  condition. 
This  is  the  volume  plasmon  radiation  composed  of  transition  radiation 
and  bremsstrah lung  radiation  that  was  described  briefly  in  the 
introduction.  It  produces  a  peak  or  shoulder  at  the  volume  plasmon 
energy  (3.70  eV  or  3300  X).  The  remaining  radiation  which 
dominates  most  of  the  spectra  in  Figs.  3  through  ",  is  roughness- 
coupled  emission  from  the  decay  of  surface  plasmons. 

The  different  shapes  of  RCSI’R  spectra  are  related  primarily  to 
the  roughness  coupling  of  plasmons  to  photons  and  not  to  the 
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production  process  of  the  plasmons.  The  most  general  qualitative 

conclusion  that  we  draw  is  that  any  mechanism  that  produces  a  rougher 

surface  has  two  effects.  It  increases  the  intensity  of  the  RCSPR  and 

shifts  its  maximum  to  longer  wavelengths  and  lower  frequencies.  In 

Fig.  3,  where  Ag  was  rapidly  evaporated  onto  a  specularly  reflecting 

Ag  substrate,  the  surface  plasmon  peak  is  relatively  weak  and  has  its 

maximum  near  the  asymptotic  energy,  trtos ,  of  the  surface  plasmon 

°  I  o  2  ) 

dispersion  curve  for  Ag  (tnus  =  3.66  eV,  >  =  3590  A).  In  Fig.  4, 
where  the  emission  is  from  Ag  on  a  matte  finish  Mo  substrate,  the 
dominant  surface-plasmon  peak  is  centered  at  a  longer  wavelength. 

We  will  argue  below  that  this  indicates  a  rougher  surface.  In 
Figs.  5  and  6,  the  very  slow  evaporation  of  Ag  results  in  increasingly 
rougher  surfaces.  Such  an  increase  in  roughness  with  film  thickness 
has  been  observed  previously. (20>24> 25)  The  effect  seems  to  be 
enhanced  in  our  experiments  when  the  electron  beam  is  incident  on 
the  radiating  region  as  new  Ag  is  deposited. 

A  satisfactory  explanation  of  the  above  results,  and  particularly 
of  Fig  s.  5  and  6,  must  account  for  the  increase  in  intensity  and 
shift  in  spectral  maximum  of  the  RCSPR  that  occurs  as  the  roughness 
is  increased.  It  should  also  account  for  the  spectral  width  and 
the  fact  that  a  nearly  constant  fraction  of  the  radiation  is 
observed  at  frequencies  above  the  asymptotic  surface  plasmon 
frequency  ws  even  after  the  spectral  maxima  has  shifted  to  a  lower 
frequency. 
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In  Fig.  8,  we  show  emission  spectra,  from  a  rough  Ag  film 
as  a  function  of  incident  electron  energies.  Smooth  surface  volume 
plasmon  radiation  is  negligible  in  the  spectra.  The  curves  are 
normalized  to  the  incident  electron  current.  In  rig.  9,  the  peak 
emission  intensity  is  plotted  as  a  function  of  electron  energy.  The 
data  points  are  fitted  by  a  theoretical  curve  discussed  in  the  following 
section. 

Fig.  10  shows  emission  spectra  for  P-  and  S-polarized  radiation 
(P  polarization  has  the  E  vector  in  the  plane  determined  by  the  inci 
dent  beam  and  the  sample  normal)  from  a  Ag  film.  The  spectra  have 
been  corrected  for  the  sensitivity  of  the  monochromator  and  optical 
system  to  the  two  polarizations.  In  the  results  shown,  the  electron 
beam  was  incident  at  5°  to  the  sample  normal,  and  the  photons  were 
detected  at  30°  from  normal.  Similar  results  were  obtained  with 
electrons  incident  at  65°  from  normal  and  the  same  observation 
angle.  This  indicates  that  the  observed  polarization  does  not 
depend  strongly  on  the  angle  of  incidence  of  the  electrons. 

These  measurements  of  the  polarization  dependence  of  the 
emission  spectra  support  our  assertion  that  two  separate  emission 
mechanisms  contribute  to  the  observed  spectra.  The  P-  and  S-polarized 
spectra  are  strikingly  different  with  the  P-polarization  curve 

O 

peaking  near  the  volume-plasmon  wavelength  (3300  A)  and  the  S-polariza- 

O 

tion  curve  peaking  near  the  surface  plasmon  wavelength  (3440  A). 

The  dashed  curve  is  the  P-S  difference  curve  and  represents  the 
excess  of  P  over  S  emission. 
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The  intensity  of  the  volume  plasmon  peak  using  smooth- surface 
transition-radiation  theory  should  increase  from  zero  the  direction 
of  the  sample  normal  to  a  maximum  about  30°  away  from  the  sample 
somal,  all  radiation  being  P-polarized.  Bremsstrahlung  is  less 
completely  polarized  since  electrons  are  randomly  scattered  in  the 
surface.  It  is  not  clear  in  the  present  data  whether  all  of  the 
excess  P  radiation  at  longer  wavelength  should  be  attril  ited  to  the 
loi  g  t :i !  1  that  is  .town  to  exist  for  the  smooth-surl ace  transition 
radiation  or  whether  some  of  it  may  represent  an  excess  P-pola.'izeJ 
emision  component  in  the  RC5PR. 

Ritchie has  derived  expressions  for  the  polarization  dependence 
cf  roughness -coupled  emission  for  plasmons  produced  by  high-energy 
grazing-incidence  electrons.  The  results  are  complex  and  stron. !v 
dependent  on  the  roughness  spectrum  but  do  indicate  an  excess  of 
P-polarized  radiation.  It  our  experiments,  electrons  are  incident 
near  normal  and  plasmons  are  produced  by  both  the  incident  electrons 
and  by  scattered  electrons  with  randomized  directions.  These  factors 
all  serve  to  reduce  the  probability  of  generating  plasmons  with 
k-vectors  in  the  plane  of  incidence  of  the  electron  beam,  which 
will,  in  turn,  reduce  the  excess  of  P  to  S  polarized  light.  For 
the  present  experiment,  where  the  electron  beam  is  incident  near 
normal,  it  appears  that  nearly  all  of  the  excess  P  radiation  can 
he  assigned  to  the  smooth  surface  volume  plasmon  radiation. 
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Qualitatively,  then,  the  difference  curve  can  probably  be  taken  as 
representing  smooth-surface  volume  plasmon  radiation  and  the 
S-polariced  curve  as  representing  RCSPR. 

Referring  to  Fig.  10,  we  see  that  the  ratio  of  the  maxima  of 
the  P-S  curves  to  the  maxima  of  the  S  curves  decreases  with 
decreasing  elect  on  energy.  In  practice  we  found  it  impossible 
to  monitor  the  smooth-surface  transition  radiation  at  electron 
energies  much  below  500  eV.  Attempts  to  increase  the  emission 
intensity  by  increasing  the  electron  flux  roughened  the  film  and 
enhanced  the  sur face-pl asmon  radiation  at  the  expense  of  the 
smooth-surface  transition  radiation. 

Discussion  of  Results 

We  wish  to  compare  the  experimental  results  on  spectral  shape  and 

(28-) 

intensity  with  the  theory  of  RCSPR  outlined  in  the  previous  paper. 

Much  of  this  discussion  is  necessarily  qualitative  since  we  do  not  have 
independent  measurements  of  the  frequency  spectrum  of  the  roughness 
auto-correlation  function.  Some  proposed  experiments  for  obtaining 
independent  measurements  of  this  function  are  discussed  briefly 
at  the  end  of  the  paper.  For  the  present,  however,  any  knowledge 
we  may  gain  about  the  roughness  spectra  must  be  extracted  from 
the  emission  spectra  themselves.  A  similar  approach  has  usually 
been  taken  by  workers  measuring  the  effect  of  roughness  in 
producing  a  dip  in  the  reflectance  of  metals  near  the  surface 
plasmon  frequency.  There,  the  decrease  in  reflectance  becomes 


larger  and  shifts  to  lower  frequencies  as  the  surface  becomes 
rougher. 

The  usual  explanation  of  the  increase  in  emission  intensity  with 
increasing  ’■oughness  is  that  the  vertical  roughness,  as  measured  by 
the  r.m.s.  deviation  of  the  true  surface  from  the  ideal  smooth 
surface,  is  increased.  This  is  the  quantity  < s  >  appearing  in  the 
theory  paper. 

The  shift  of  the  emission  peak  to  lower  frequencies  (longer 
wavelengths)  might  be  explained  by  attributing  the  radiation  to 
relatively  small  q  roughness  components  which  couple  plasmons  in 
the  retarding  region  to  the  radiation  field  and  by  assuming  that  the 
roughness  spectrum  shifts  to  lower  q  values  with  increasing  roughness. 

But  this  mechanism  is  not  consistent  with  our  observation  of  strong 
radiation  above  <vs ,  nor  does  such  an  assumption  lead  to  reasonable 
predictions  of  the  dependence  of  the  radiation  intensity  on  electron 
energy.  In  the  retarding  region,  (k'2)  equation  (2)  shows  that 
plasmon  production  is  (.nearly)  linear  in  electron  energy  in  contrast  ro 
the  data  shown  in  fig.  9.  The 

second  possibility  for  explaining  the  shift  to  lower  frequencies 

as  t!  it  the  peak  radiation  is  produced  by  high  q  roughness  components, 

but  that  increasing  roughness  somehow  decreases  the  frequency  of  high 

k- vector  plasmons.  Other  workers  have  also  obtained  experimental 

results  which  they  attribute  to  a  shift  of  the  surface  plasmon 

curve  to  lower  frequencies  (or  equivalently  to  a  reduction  in 

plasmon  phase  velocities)  as  the  surface  is  roughened .  (2 1 )  An  alternative 

'‘V ! an. at  ion  for  this  shift  can  be  given  by  assuming  that  the  rough  surface 
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acts  as  a  surface  layer  with  a  refractive  index  intermediate  between  that  of 
air  and  Ag,  but  to  our  knowledge  no  rigorous  theory  accounting  for  the  shift 
has  been  given. 

The  attractive  feature  of  the  calculations  of  Kretschinann  e^  ad .  is  that 
they  provide  explanations  for  the  major  features  of  the  RCSPR  spectra 
both  as  a  function  of  electron  energy  and  as  a  function  of  increasing 
vertical  roughness  using  standard  surface  plasmon  and  surface  roughness 
theory  and  the  single  additional  assumption  that  the  roughness  spectrum 
has  a  prominent  peak  at  high  q-vector  that  is  responsible  for  the 
RCSPR  in  the  vicinity  of  its  spectral  peak. 

For  simplicity,  the  high  q-vector  roughness  peak  has  been  modeled 
as  a  delta  function.  Fig.  la  of  the  theory  paper  shows  the  decrease 
of  the  ratio  of  the  RCSPR  to  the  volume  plasmon  emission  with  increasing 
electron  energy.  This  is  in  agreement  with  our  observation  that 
RCSPR  increasingly  dominates  the  radiation  at  low  electron  energies. 

Fig.  2  of  the  theory  paper  shows  the  emission  spectrum  as  a 
function  of  q2  <s?>  where  is  the  dominant  roughness  component  and 

<s2>  the  r.m.s.  height  of  the  roughness.  The  surprising 
result  is  that  the  single  high-q  roughness  component  produces  a 
double  peaked  radiation  spectrum  with  one  peak  lying  above  us  and 
the  second  and  much  larger  peak  lying  below  ws .  Both  the  intensity 
and  the  splitting  increase  with  increasing  <s2>.  If  the  delta 
function  in  the  roughness  spectrum  is  replaced  by  a  high  q-vector 
peak  of  finite  width,  additional  broadening  is  introduced  so  that 
the  splitting  is  no  longer  so  well  resolved.  The  shift  of  the 
dominant  peak  position  to  lower  frequencies  and  larger  wavelengths 
with  increasing  ss2--  is  retained  however.  Similar  splittings 
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but 


have  been  noted  for  experiments  on  diffraction  gratings, 
the  effects  are  much  larger  for  large  q-vector  roughness  peaks 
since  the  splitting  scales  with  a2. 

The  splitting  is  not  resolved  in  Ag,  which  is  not  surprising 
since  the  model  does  not  include  a  finite  width  for  the  roughness 
spectrum  nor  does  it  treat  the  lifetime  broadening  of  the  surface  plasmon  in 
Ag  adequately.  It  would  be  difficult  to  observe  a  splitting  for 
Ag  in  any  case,  because  u>v  and  ws  occur  very  close  together.  But 
the  theory  does  provide  suggestive  explanations  for  the  dependence 
of  peak  radiation  intensity  on  electron  energy,  and  for  the  shift 
in  peak  position  as  roughness  is  increased.  Neither  of  these 
results  of  the  calculations  is  changed  significantly  if  a  roughness 
peak  of  finite  width  is  assumed. 

As  indicated  in  Eq.  6  and  Fig.  lb  of  Kretschmann  et  al. 
the  peak  intensity  versus  electron  energy  is  scaled  along 
the  horizontal  axis  by  q^':  where  q^2  designates  the  dominant 
roughness  component.  Using  this  universal  curve,  we  can  fit  the 
experimental  values  of  the  peak  intensity  versus  electron  energy 
shown  in  Fig.  9  if  we  assume  that  q^  =  1.12  x  10~?  A  !,  which 

O 

corresponds  to  a  wavelength  of  =  560  A,  for  the  dominant  roughness 
component.  This  is  equivalent  to  a  value  of  q^  =  6.25  in  the  reduced 
units  of  (w/c)  that  were  used  in  the  introduction  of  this  paper  and 
by  Kretschmann  et  al .  As  we  saw  in  the  introduction  the  dependence  of 
peak  intensity  on  the  electron  energy  results  primarily  from  the  energy 
dependence  of  t he  probability  of  producing  plasmons  of  wavevector  kyq^. 
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Fitting  the  data  of  Fig.  9  with  Fq.  2  of  this  paper  gives  very 
nearly  the  same  value  of  (q^  ^  k  =  8)  as  the  more  exact  treatment. 

The  Kretschmann  £t_  theory  also  predicts  a  shift  in  the  dominant 
peak  position  which  is  analytically  reiated  to  the  change  in  peak 
height.  As  may  be  seen  from  Eqs.  6  and  7  of  the  theory  paper,  a 


plot  of  emission  intensity  versus  4[(1 


op)/(l  '  cp)]‘  >  where 


is  the  location  of  the  spectral  peak  measured  in  terms  of  the  real 


part  of  the  dielectric  constant,  should  give  a  straight  line  as  the 
r.m.s.  height  <s2>  of  the  dominant  roughness  component  is  increased. 
This  has  been  done  for  the  data  of  Fig.  6  and  is  shown  in  Fig.  11. 


The  observed  straight  line  plot  provides  support  for  our  contention 
that  a  single  dominent  high  q-vector  component  of  roughness  can 
account  satisfactorily  for  most  of  our  data. 


There  are  some  independent  indications  that  a  large  high  q-vector 
peak  occurs  in  the  roughness  autocorrel at  ion  function.  The  growth 
process  can  easily  lead  to  such  a  result,  since  these  evaporated 
films  usually  develop  from  isolated  islands  which  grow  with  the 
addition  of  material  until  they  touch .  (2 1 > 34)  It  scons  reasonable 
that  there  is  a  dominant  horizontal  roughness  component  whose  scale 
is  set  by  the  average  size  of  the  islands  when  they  first  touch, 
and  in  fact  most  photomicrographs  of  sufficient  resolution  show 
regular  nucrostructure  with  lateral  dimensions  of  a  few  hundred 
A. (  ■ • "  )  Autocorrelation  functions  obtained  from  interferometric'  1 
and  light  scattering*^'  ^  neasurenent s  .ire  sensitive  primarily  to  small 
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q  roughness  components  and  can  not  give  evidence  of  a  high-q  roughness 
peak.  The  photomicrographic  technique  of  Dobberstein  does  indicate 
the  presence  of  a  high-q  peak  in  the  roughness  spectrum.  For  Ag , 
Dobberstein  found  a  peak  at  7  x  10"  A’1  and  kaspar  and  Kriebigt-’1) 
found  peaks  at  7  x  10  A  1  and  15x10  A  1 .  Braundmeier  and 
Hall  found  it  necessary  to  assume  an  autocorrelation  function  with 
a  peak  at  2  x  10  A  1  in  order  to  explain  the  radiation  patterns 
of  optically  excited  surface  plasmons  on  Ag  films. These  values 

O 

may  be  compared  with  the  value  of  11  x  10  A  obtained  from  fitting 
the  data  of  Fig.  9. 

It  would  be  desirable  to  perform  experiments  on  surfaces  with 
completely  characterized  roughness  functions.  By  combining  methods 
sensitive  to  small  q^-5-  ’ components  with  photomicrographic  methods^'  ^ 
sensitive  to  high-q  components,  it  may  be  possible  to  specify  the 
roughness  spectrum  over  the  full  region  of  interest.  To 
confirm  the  splitting  predicted  theoretically,  it  would  be  better 
still  to  prepare  surfaces  with  a  single  dominant  q-vector  as  has 
been  done  for  small  q-vectors  using  diffraction  gratings. We 
hope  to  pursue  both  approaches  in  the  future. 
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FIGURE  CAPTIONS 


Fig .  1 . 

Fig.  2 . 
Fig.  3. 


Fig.  4. 
Fig.  5. 

Fig.  6. 

Fig.  7. 

Fig.  8. 
Fig.  9. 

Fig.  10 


(a)  Dispersion  curve  for  surface  plasmons 

(b)  U'avevector  conservation  for  roughness-coupled  surface- 
plasmon  radiation 

Schematic  of  experimental  apparatus 

Emission  spectrum  for  Ag  on  Ag  showing  both  surface-  and 
volume- pi asmon  peaks  ( —  before  evaporation;  —  after 
evaporation.)  •  Volume-plasmon  and  limiting  surface -pi asmon 
wavelengths  are  indicated  on  this  and  most  subsequent 
figures . 

Emission  spectrum  for  Ag  on  matte  Mo  substrate 

Growth  of  low-energy  emission  peak  with  slow  eva. -oration 

of  additional  Ag  on  Ag 

Growth  of  emission  intensity  with  Ag  thickness  for  deposits 
of  Ag  on  glass 

Shift  of  emission  peak  as  the  roughness  of  a  Ag  film  i 
increased  by  electron  beam  heating 

Emission  spectra  for  various  inciden*  electron  energies 
Peak  intensity  as  a  function  of  electron  energy.  Smooth 
curve  is  from  theory  of  Kretschman  et  a_l  described  in  text. 
P  and  S  polarized  emission  spectra  from  relatively  smooth 
surface  showing  both  volume-plasmon  radiation  and  roughness- 
coupled  surface- p lasmon  radiation. 
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Fig  11.  Plot  of  peak  intensity  versus  4 [ ( 1  +  e)/(l  -  r)]  ,  a 

function  of  peak  position.  The  plot  should  be  a  straight 
line  in  the  theorv  of  Kretschnann  et  al. 
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CHAPTER  II-e 


Theory  of  Emission  Spectra  from  Metal  Films  Irradiated  by  lo w  Energy 
Electrons  near  Normal  Incidence* 

"f"  "f  ^ 

E.  Kretschmann,  T.  A.  Callcott,  ‘  and  1.  T.  Arakawa 
Health  and  Safety  Research  Division,  Oak  Ridge  National  Laboratory 
Oak  Ridge,  Tennessee  37830 

ABSTRACT 

The  emission  spectrum  produced  by  low  energy  electrons 
incident  on  a  rough  metal  surface  has  been  calculated  for  a 
roughness  auto-correlation  function  containing  a  prominent  peak 
at  a  high  wavevector.  For  low  energy  electrons  near  normal 
incidence,  the  high  wavevector  peak  dominates  the  roughness 
coupled  surface  plasmon  radiation  (RCSPR)  process  The  cal 
cnlation  yields  estimates  of  the  ratio  of  RCSPR  to  transition 
radiation,  the  dependence  of  emission  intensity  on  electron 
energy  and  the  shape  and  position  of  the  RCSPR  peak.  The 
most  interesting  result  is  that  the  high  wavevector  roughness  car.  split 
the  RCSPR  radiation  into  peaks  lying  above  and  below  the 
asymptotic  surface  plasma  frequency.  The  results  are  compared 
to  data  from  Ag  in  the  following  paper. 
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I.  INTRODUCTION 


A  large  number  of  experimental  and  theoretical  studies  of 
the  emission  of  light  from  metal  surfaces  irradiated  by  elections 
have  been  performed  in  the  last  two  decades.  The  most  recent 
ones  are  cited  here;  the  earlier  ones  are  listed  in  references 
1-4.  In  the  emission  spectra  an  edge  or  weak  maximum  is 
generally  found  experimentally  neai  id  =  (.the  volume  plasma 
frequency  of  the  metal).  Under  certain  conditions  an  additional 
laximum  is  found  near  u  =  (the  surface  plasma  frequency). 
Theoretically  'he  edge  at  u.^,  can  be  described  by  transition 
radiation  (and/or  Bremsstrahiungj  at  a  smooth  surface/’  and  the 
maximum  near  by  roughness  coupling  of  nonradiative  "surface 
plasmons"  to  the  radiation  field/’ 

Ther:  are  three  published  theoretical  treat¬ 
ments  of  the  latter  effect.  In  one  quantum  mechanical  calculation 

4 

and  in  a  classical  calculation,  based  on  earlier  theoretic  ■] 

8 

work,  the  electrons  are  assumed  tc  have  nearly  gracing  incidence 
and  only  the  radiation  from  high  wavevector  surface  plasmons  are 
calculated.  The  change  of  the  dispersion  relation  of  surface 
plasmons  due  to  the  roughness  is  not  included.  A  more  general 
calculation  has  been  carried  out  recently'^  which  involves  all 
angles  of  incidence  as  well  as  the  influence  of  surface  roughness 
on  the  dispersion  relation.  But  in  that  study  only  the  case 
of  nearly  gracing  incidence  electrons  was  discussed  in  detail 
because  till  now  most  experiments  have  been  done  in  this 


v  ■ 
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configuration  in  order  to  enhance  the  surface  plasmon  radiation. 

Recently  transition  radiation  experiments  at  rough  metal  surfaces 
have  been  performed  with  low  energy  electrons1  (E^  SI  keV) .  It 

was  found  that  the  natural  roughness  of  evaporated  metal  films 
causes  the  intensity  of  the  surface  plasmon  radiation  to  be  much 
great°r  than  that  of  the  transition  radiation  even  for  normally 
incident  electrons.  We  include  in  the  calculation  the  effect  of 
the  roughness  on  the  dispersion  relation,  which  has  been  shown  to  cause 
splitting  of  the  high  wuvevector  part  of  the  dispersion  relation  into  two 
branches . ^ 

The  model  used,  the  results  of  the  calculations,  and  a 
discussion  of  the  results  for  a  particular  correlation  function 
will  be  given  here.  A  comparison  of  experimental  data  on  Ag  with 
this  theory  will  be  given  in  the  following  paper. 

II.  MODEL  AND  RESULTS 

An  electron  of  charge  -e  penetrates  the  metal -vacuum 
boundary  at  an  angle  a  with  the  surface  normal  with  a  constant 
velocity  v.  Electron  velocity  is  small  compared  to  the  velocity  of  light 
(B  =  v/c  <<  1) .  The  metal  is  described  by  the  dielectric 

function  e(w).  For  a  smooth  surface  the  emission  spectrum  is  calculated  as 
follows.  The  homogeneous  fields  associated  with  the  moving  electron  in 
both  vacuum  and  the  metal  are  calculated  by  the  Maxwell  equations. 
Homogeneous  fields  in  the  metal  and  the  vacuum  are  added  in 
order  to  satisfy  the  usual  boundary  conditions  at  an  interface. 

The  fields  have  radiative  parts  (transition  radiation)  where  the 
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wavevector  along  the  boundary  is  smaller  than  w/c,  the  wavevector 
of  light,  and  nonradiative  parts  (surface  plasmons).  Additional 
details  of  the  method  and  the  results  of  the  smooth  surface 
calculation  can  be  found  in  ref.  10.  From  these  zero  order  fields, 
which  are  present  at  a  smooth  surface,  the  scattered  fields 
due  to  the  surface  roughness  up  to  the  second  order  in  the 
roughness  height  are  calculated  with  a  method  given  in  ref.  11. 

From  the  zero  order  fields  we  get  for  the  power  of  the 
emitted  electromagnetic  radiation  in  the  plane  of  incidence  per 
frequency  interval,  dw,  observed  at  the  angle  9  within  an  element 
of  solid  angle  dli 

iLJ —  =  ( — % — )  cos2  a  sin2  9  cos2  e ) e - 1 j  2  02/ ]e  cos  9  +  /c-sin29  |2. 
dfiduj  irzc  i 


The  Hq.  (1)  can  be  found  in  ref.  5  also.  It  predicts  a  weak 
maximum  near  the  plasma  frequency  when  e  =  0,  which  is  usually 
called  the  "volume  plasmon  peak."  From  the  scattered  fields  we 
get  for  the  special  case  of  normally  incident  electrons  (a  =  0) 
and  for  the  detection  angle  9=0 


d2!1  ,  ,  e 


dfiduj  ~  ^  m':C  ' 

fdk  k3  gf-4- 


-  1  |  2  (— )4  <s2> 

I  C 


k  i  +  k  . 


k> 


S2 

(i  +  k2e2)2 


kj  ♦  ek2  1 

Mere  kj  =  -  k2  ,  k2  =  I  k2  (Imki ,  Imk2  >  0)  and 

k  is  the  wavevector  of  the  zero  order  waves  along  the  surface 
in  units  of  u>/c.  The  surface  properties  enter  through  <s2>, 


(2) 
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the  mean  square  roughness  height,  and  g(— •  k)  the  normalized 
correlation  function  of  the  roughness.  (For  a  definition  of 
g(— •  k)  see  ref.  11. J  In  Eq.  (2)  terms  with  02  (but  not  with 
S2k2)  are  neglected  and  only  the  transverse  magnetic  (TM)  zero 
order  fields  are  taken  into  account.  The  TM  fields  are  associated 
with  surface  plasmons  and  are  of  much  larger  amplitude  than  transverse 
electric  (TE)  fields. 


III.  DISCUSSION 

We  will  discuss  Fq.  (2)  for  frequencies  where  surface 

plasmons  may  be  present  (e(u)  0.)  For  the  discussion  we  will 

use  a  special  roughness  correlation  function  with  features  that 

we  believe  roughly  approximate  the  true  correlation  function. 

The  correlation  function  of  a  rough  ^Tirface  can  be  determined 

in  this  low  roughness  k-vector  region  by  light  scattering 
12 

experiments.  It  has  been  generally  assumed  that  the  correlation 
function  has  a  strong  maximum  for  k  =  0  and  falls  monotonically 
with  increasing  k(w/c).  However,  there  are  indications  that  in 
the  high  k-region  the  correlation  function  has  a  second  maximum 
at  some  k  -  k^  >>  1  Here  k^  =  ig/X^  with  X^  the 

dominant  wavelength  of  the  statistical  roughness  and  X^  the 
wavelength  of  light  at  the  surface  plasma  frequency  a >  .  Therefore, 
we  divide  g(— •  k)  into  a  low  k-vector  part,  g^(kl,and  a  high 
k-vector  part,  g^(k)  which  will  be  approximated  by  a  6-function. 
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<s2>  g(—  •  kj  =<s  2>  g.  (,~  •  k)  +  ^r-  <s,2> - - -  5  (.k  -  k.  )  ,  (3] 

c  *.  *■  c  m  n  n 

V  Kh 

where  <s2>  =  <s  2>  +  <s, 2>  and  k,  =  A  /A.  .  More  generally  we  mav 
l  h  h  s  h  & 

introduce  a  sum  of  6-functions  to  describe  the  high  k-vector  part. 

k2  +  k  j  |  2 

Equation  (2)  can  be  simplified  since  — r - —  |  has  a  pole 

C  K2  +  K \  I 

if  k  =  k  =  /  e/(e  +  1)  ,  the  normalized  k-vector  of  surface 
sp 

plasmons.  For  the  scattered  power  of  "low  k-vector"  surface 
plasmons  (e  <  -l.S,  k  £2,  Ijtie<<  Re  e) ,  we  approximate  the 
integral  in  Eq.  (2)  by  the  integral  over  the  pole  and  obtain 


d2I  ,  7  ,  e2  ,  1  -  Re  e 

— —  =  4irz  f — 7 - j  •  — - - 

duidft  '  ir'c  ime 


k0G  g(-f-  •  k0)  14) 


where  kQ  =  /  Re  e/(.Re  e  +  1J  .  Equation  (4)  is  similar  to  the 
equation  which  has  been  developed  to  describe  the  dip  in  the 
reflection  near  the  surface  plasma  frequency.  For  the  scattered 
power  of  "high  k-vector"  surface  plasmons  (e  =  -1,  k  »  1)  we  get 

-£sr  “ 8*  <sh2>  • 


/e  -  ll2 


E  ♦  11 


dk  k3  gh(-f- 


(MkW 


In  order  to  get  Eq.  (.5)  we  set  ik  =  k  x  =k2  which  is  appropriate 
for  large  k  With  the  special  choice  of  g(-^—  •  k  )  in  Eq.  (.3)  we  get 


1 
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aZ  *  <Sh2>  V 


1 

8.2 


(8) 


For  further  discussion  we  will  perform  some  numerical  calculations 
which  treat  aspects  of  the  theory  that  can  be  compared  to 
experiment.  We  calculate: 

1.  The  intensity  relation  between  "surface  plasmon" 
and  "volume  plasmon"  peaks. 

2.  The  electron  energy  dependence  of  the  surface 


1 . 


plasmon  peak  intensity. 

3.  The  shape  (half  width  and  maximum  position)  of 
the  surface  plasmon  peak. 

4.  The  relation  between  maximum  position  and  intensity 
of  the  surface  plasmon  peak. 

For  0  =  30°,  a  =  0,  e  =  0,  Eq .  (1)  gives  J°  =  = 

Cl  1  id  (l) 

3  2 

—  -7 —  B2  as  the  maximum  of  the  volume  plasmon  peak.  The 

4  TT  c 

ratio  of  the  intensity  of  the  surface  plasmon  peak  (  Eq .  (6) 
for  e  =  -1  +  ie2)  and  the  volume  plasmon  peak,  is  given 
in  Fig.  la.  Figure  la  shows  clearly  that  a  small  B  value 
(or  a  low  electron  energy)  is  necessary  in  order  to  enhance 


the  surface  plasmon  peak  relative  to  the  volume  plasmon  peak. 


If  we  take  E  =  1  keV  and  Ime  =  e2  =  0.3  we  need 
<s2>l/2/^  >  1/J52  to  get  J  >  10  J°.  If  we  took  E  =  10  keV 

we  would  get  J^1  %  .1°  for  the  same  roughness  parameter. 
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The  energy  dependence  of  the  surface  plasmcn  peak  intensity  is 

shown  in  Fig.  lb.  For  small  is  a  linear  function  of  B7 .  It 

peaks  for  £  =  1/k,  .  If  one  fits  the  curvature  of  the 
h 

experimentally  determined  energy  dependence  with  the 
theoretical  one,  one  can  calculate  the  dominant  wavelength 
of  the  roughness. 

3.  In  Fig.  2  roughness  coupled  surface  plasmon  spectra  as  predicted  bv 

F.q.  (61  with  the  replacement  (7)  are  shown.  The  Ret  values  of  silvc 

determined  in  ref.  13  are  assumed  and  Ime  is  set  equal  to 
1 8 

0.3.  Spectra  for  different  values  of  the  roughness 

parameter  b  =  /<s^>  /\  are  plotted.  The  denominator 

(1  +  k  •' Q 2 ) ?  has  been  neglected  Whereas  one  peak  is  predicted  near 

e  =-l  for  b  =  0,  a  double  peak  is  predicted  for  b  '  0.  The 

high  energy  peak  is  always  smaller  in  intensity.  The 

splitting  increases  with  b.  Recently  this  double  peak 

structure  was  found  '  ^  in  the  emission  spectrum  of  potassium. 

an 

There  is  also  a  report  in  the  literature"  which  shows 
double  peaks  in  reflection  measurements  of  sodium  and 
potassium.  Because  the  splitting  changes  with  the  value  of 
k  it  may  he  smeared  out  since  in  real  physical  systems 
g^  [  (ui/c)k]  contains  more  than  one  Fourier  component.  A 
broad  peak  which  is  shifted  to  smaller  e-value  (larger 
wavelength)  may  be  seen  in  this  case.  This  is  the  situation 
that  we  believe  to  he  appropriate  to  the  Ag  data  presented 
in  the  following  paper. 
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4.  In  Fig.  2  it  can  be  seen  that  an  increase  of  the  intensity 

is  always  connected  with  a  shift  in  the  position  of  the  maximum 

to  longer  wavelength.  According  to  Eq.  (6)  the  intensity 

is  proportional  to  anu  according  to  Eqs .  (7)  ana 

(8)  ie;  -  l)/(i.  +  1)  =  a  is  proportional  to 

/<s  2>k.  2  .  This  result  should  be  vajid  for  all  experi- 

h  h 

mental  situations  where  roughness  coupling  to  high  k-vector 

surface  plasmons  takes  place.  This  relation  between 

intensity  and  shift  vas  found  in  the  reflection  measurement 
2 1 

at  rough  surfaces,  and  is  demonstrated  for  emission  from 
Ag  ii  the  following  paper. 

The  results  obtained  in  the  numerical  calculations 
associated  with  paragraphs  1  through  4-, above  are  consistent  with 
experimental  results  on  emission  spectra  from  rough  silver  films 
reported  in  the  accompanying  paper. 
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FIGURE  CAPTIONS 


Fig.  la  Ratio  of  the  intensities  of  the  "surface  plasma" 
peak  to  the  "volume  plasma"  peak  Ju. 

Fig.  lb  Intensity  of  the  "surface  plasma"  peak  J^1  versus 
incident  electron  energy. 

Fig.  2  Wavelength  dependence  of  the  emission  intensity  J 
(Eq.  6,  7,  8)  from  high  k-vector  surface  plasmons 
The  values  of  Re  £  are  shown  below  the  wavelength 
scale;  Ime  is  assumed  to  be  0.3.  The  curves  are 


labeled  by  values  of  the  relevant  roughness  parameter 


74 


CHAPTER  II-C 


A  SPLITTING  OF  THE  DISPERSION  RELATION 
OF  SURFACE  PLASMONS  ON  A  ROUGH  SURFACE* 

+  +t 

E.  Kretschmann,  T.  L.  Ferrell,  '  and  J.  C.  Ashley 
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Abstract 

The  dispersion  relation  for  surface  plasmons  on  a 
statistically  slightly  rough  surface  is  shown  to  display 
a  splitting  under  certain  conditions.  For  a  simple  model 
of  the  correlation  function  of  the  surface  roughness, 
explicit  calculations  are  made  which  explain  experimental 
results  not  previously  interpreted. 
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The  dispersion  relation  for  surface  plasmons  on  a  statistically 
slightly  rough  surface  is  of  continuing  interest  in  surface  physics  as 
attempts  are  made  to  characterize  more  closely  actual  surfaces.  The 
present  paper  considers  this  dispersion  relation  for  surface  plasmon 
wave  vectors  K  large  relative  to  uj/c,  where  w  is  the  surface  plasmon 
frequency. 

The  dispersion  relation  for  surface  plasmons  on  a  smooth  solid-vacuum 
boundary  can  be  written 

ns(K)  =  £  <2  -f  K1  =  e  Vk2  -  1  +  Vi'2  -  c  =  0,  (la) 

or  simply 

1_ 

c  +  1=0.  (lb) 


In  these  equations,  K  is  the  (magnitude  of  the)  surface  plasmon  wave  vector 
in  units  of  w/c.  The  function  e  =  e (c^ )  is  the  complex  dielectric  function 
of  the  solid.  The  dispersion  relation  is  shown  in  Fig.  1  for  In  c  =  0. 

If  the  surface  is  rough,  the  dispersion  relation,  Eq.  (1),  must  be 

modified.  Equations  for  the  modified  dispersion  relation  have  been  given 

1  -4 

by  several  authors.  The  most  convenient  form  for  the  present  purposes 

3 

is  given  by  Toigo,  Marvin,  CeT 11,  and  Hill.'  We  have  verified  the 
equivalence  of  their  results  to  those  of  Kroger  and  KretschmannJ 
The  dispersion  relation  for  a  surface  with  statistically  small  roughness 
may  be  written  as 

2 

nR(K)  =  n$(K)  -  -  l)?<s2>  I(K)  =  0,  (2) 
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,  * 

where  /  is  the  "mean- square  roughness  height,"  and 


I(K)  =  f  d2K'g(^|  K  -  K'|) 

r  J 


(K-iT  -  K2Kf)  (&•&' 

lOT'l 


-  k2'K1  ) 


(3) 


with  g(~jK  -  K'|)  the  surface-roughness  correlation  function  K-f  = 

K-|  (K') ,  and  =  KglK') . 

Equation  (2),  or  an  equivalent  form,  has  been  discussed  quali¬ 
tatively  in  connection  with  the  broadening  of  the  surface  plasmon 

resonance  (increased  damping),  and  the  increase  of  the  wave  vector  of 

1  c 

"low  wave  vector"  surface  plasmons.  *  The  "low  wave  vector"  region 
corresponds  tc  K  '  2.  Good  agreement  of  these  expectations  with  the 

g 

experimental  changes  in  the  dispersion  on  a  grating  has  been  found. 

In  the  "high  wave  vector"  region,  K  >  2,  the  dispersion  relation  for 
a  smooth  surface  simplifies  to  e(u)  as  -1,  as  can  be  seen  from  Eq.  (lb) 
or  Fig.  1.  When  surface  roughness  is  included,  a  new  feature  arises  in 
the  dispersion  relation  in  this  wave  vector  region,  in  addition  to  the 
obviously  different  dependence  on  K,  which  prevents  a  simple  extrapo¬ 
lation  of  tne  results  described  for  K  <  2  to  values  of  K  >  2.  This  new 
feature  in  the  "high  wave  vector"  region  is  a  splitting  of  the  dispersion 
relation  for  surface  plasmons  on  a  rough  surface.  This  effect  is  demon¬ 
strated  and  described  in  the  remainder  of  this  paper. 

If  the  surface  wave  vector  has  a  large  enough  magnitude  so  that 
one  neglects  terms  of  order  K  compared  to  K,  one  may  put  as  x  K 
in  F.q.  (la)  to  obtain 

ns(K)  »  K(e  +  1).  (4) 
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For  the  present  purposes  we  divide  I(K)  of  Eq.  (3)  into  two  integrals, 
1^  and  12,  given  by  I  =  K [i ^ / ( e  +  1)  +  1^] ,  where  (1^  is  restricted 
the  region  K'  <  2  (K'  >2).  In  the  "high  wave  vector"  region  of  I,  n 
is  approximated  by  the  limiting  form  in  Eq.  (4),  and  Eq.  (2)  becomes 

r  =  * + 1  -  y- -  '^^(rrKr  >2] -  <*> 

where,  for  K  >  2, 

^  \  I  d2K'g(||K  -  K'|)KK'(1  -  cos*)2  ,  (6) 

c  * 

K'>2 

2  r  ?  (K'cos;  -  K4) (K'cosf  -  K{) 

I2*T  J  dVg(^K-  K'|)K-  - —  *  (7) 

c  K'<2  s 


1  a 
1  +  a 


(11) 


>  I,  .  (12) 
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and  g  =  1/a.  Equation  (1(J)  yields  the  splitting.  This  surprisingly 
simple  result  has  tne  consequence  that,  under  the  conditions  dictated  by 
tne  aoove  approximations,  one  may  observe  a  double  peak  near  the  surface 


plasma  resonance.  The  separation  of  the  peaks,  or  the  splitting,  is 

determined  by  the  difference  in  the  roots  of  Eq.  (10)  and  is  given  by 

- 

A e  =  p  -  a  ~  4* /  <; s‘  b  I,  =  4a. 


In  Fig.  2  we  plot  the  response  fjnction  from  Eq.  (5)  with 

I;  -  0,  as  a  function  of  Re  e .  To  be  more  specific  we  have  chosen  Im  r  = 
0.1,  corresponding  to  potassium,  and  the  wavelength  scale  at  the  top  of 

7 

the  figure  is  for  potassium  with  fuc  '  2.8  eV.  These  curves  give  the 

values  of  c(^)  at  which  surface  plasmon  excitation  is  possible  (positions 

of  the  maxima)  and  the  damping  of  these  modes  (half-width).  The  curve 

labeled  a“  =  0  corresponds  to  surface  waves  on  a  smooth  surface.  The 

curves  for  z~  >  0  show  the  splitting  found  in  the  response  function  (o^ 

the  dispersion  relation)  when  surface  roughness  is  included  in  the  descrip- 

2 

fciori  of  nigh  wave  vector"  surface  waves.  For  a  given  value  of  a‘  >  0,  the 
surface  waves  at  a  lower  e  value  (larger  vacuum  wavelength)  always  have  a 
larger  intensity  than  those  with  a  higner  e  value  (smaller  vacuum  wave¬ 
length 

In  order  to  analyze  Eqs.  (3)-{l?)  ir,  more  detail,  we  choose  a 

special  form  of  the  correlation  function  g(wK/c).  It  is  likely  that 

most  metal  samples  have  a  roughness  wavelength  which  peaks  around  a 
8  9 

given  value  X^.  *  Therefore,  we  assume 


2  *  KpW 


2  «(K  -  Kr), 


g(wK/c)  a 


(14) 


where  the  factor  in  the  denominator  of  Eq.  (14)  has  been  chosen  so 

that 

2  r 

^  I  d  K  g(a»K/c)  =  1.  (15) 

c  J 

In  general  the  correlation  function  may  be  given  by  a  linear  superposition 
of  functions  like  that  of  Eq.  (14),  and  there  is  always  a  maximum  in 
g(wK/c)  for  K  =  0. 

As  the  case  KR  <<  K  below  shows,  the  choice  of  a  form  for  g  under 
the  above  simple  conditions  has  very  little  effect  on  the  dispersion 
relation  for  long  roughness  wavelengths. 

Inserting  Eq.  (14)  into  Eq.  (6)  for  I,,  we  calculate  the  integral 
for  the  cases  KR  <<  K,  KR  =  K,  and  KR  >>  K.  These  cases  respectively 
imply  a  surface  roughness  wavelength,  AR,  larger,  equal  to,  or  smaller 
than  the  surface  plasmon  wavelength.  The  results  are  given  in  Table  1. 


Table  1.  Numerical  Examples  (A  =  4000  &, 
K  =  10,  <s2>h  =  Ar/20) 


*1 

kr 

<s2>  "7ar 

Ac 

3  4 

Kr  «  K 

2 

100/2000 

0.077 

32K2 

KR  =  K 

(1  -  2k2 

10 

20/400 

0.44 

kr  »  K 

1.5  KKr  -  2K2 

40 

5/100 

0.63 

The  table  also  contains  the  splitting  Ac  for  KR  =  2,  10,  and  40  if 
electromagnetic  waves  with  a  vacuum  wavelength  A  -  4000  are  involved 
in  the  surface  plasmon  excitation  process.  We  have  taken  the  surface 

2  ^ 

plasmons  to  have  K  =  10  and  the  surface  to  have  a  roughness  height  <s  > 
given  by  <s<>)>'5  =  AR/20. 
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The  numerical  examples  show  that  even  a  small  roughness  height  can 
produce  a  remarkable  splitting  of  the  dispersion  relation.  If  one 
wishes  to  describe  properly  the  behavior  of  electromagnetic  waves  at  a 
rough  surface  for  frequencies  near  the  surface  plasma  resonance,  the 
change  in  the  dispersion  relation  should  be  included.  However,  this  has 
not  been  done  in  the  analysis  of  reflection  measurements  on  rough  sur¬ 
faces.10’11 

The  splitting  of  the  dispersion  relation  should  be  experimentally 
observable.  There  is,  in  fact,  a  report  of  a  double  peak  near  the 
surface  plasma  frequency  in  reflection  measurements  on  sodium  and 
potassium.^  For  their  potassium  sample  this  splitting  was  observed  to 

be  approximately  0.2  eV ,  implying  that  Ac  rs  0.3.  From  our  calculations 

2  U 

above  for  K  =  Kp,  this  would  mean  <s  >2/Ap  ~  1/30,  which  is  quite 
reasonable. 

Recently,  in  the  light-emission  spectra  from  films  bombarded  by 

low-energy  electrons  1  keV),  a  doublet  structure  at  the  surface 

1 2 

plasma  energy  was  found  for  some  potassium  films.  This  structure  is 

readily  explained  in  terms  of  the  theory  presented  here. 

In  an  experiment  in  which  surface  plasma  waves  of  a  variety  of  wave 

vectors  are  involved,  the  splitting  may  not  be  detectable.  This  is  due 

to  the  fact  that  the  superposition  of  the  associated  response  functions 
2 

| K/n R i  ,  each  with  its  own  splitting,  may  smear  the  results  into  a  broad 
maximum  with  a  lower  e  "alue  (<  -1).  Additionally,  there  will  be  a 
substantial  amount  of  response  for  -1  <  e  <  0,  which  comes  from  the 
branch  at  the  higher  e  values. 
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Figure  Captions 


Fig.  1.  Dispersion  relation  of  surface  plasma  waves  on  a  smooth  surface. 

?  (  _  i\2«2r2 

Fig.  2.  The  response  function  |K/nR(  =  (e  +  1)  -  '£  ~  +  —  for 

surface  plasmons  on  a  rough  surface  (a  >  0).  The  response  function 

2 

for  a  smooth  surface  is  given  by  the  curve  labeled  a  =0.  Each 

2 

curve  is  scaled  so  that  the  maximum  value  of  |K/n^|  is  equal  to 
one. 
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Abstract 


Emission  spectra  from  potassium  films  irradiated  with  1.5-keV 
electrons  have  sfown  a  splitting  of  the  roughness-coupled  surface- 
plasmon  radiation  into  two  peaks,  as  predicted  by  a  recently  developed 
theory. 


It  has  been  shown  theoretically  [1,2]  that  high-wavevector  compo¬ 
nents  in  the  roughness  spectrum  for  a  surface  should  cause  the  roughness- 
coupled  surface-plasmon  radiation  curve  to  split  into  two  peaks  lying 
above  and  below  the  normally  observed  asymptotic  surface-plasmon  energy 
This  splitting  was  looked  for,  but  not  resolved,  in  emission  spectra 
from  election-irradiated  Ag  films  [3].  This  letter  demonstrates  its 
existence  in  the  emission  spectra  from  electron-i rradiated  "smooth"  K 
films. 

high-wavevector  components  'in  the  surface-roughness  spectrum  arise 
from  small  scale  irregularities  or  structures  which  may  be  superimpose! 
on,  and  independent  of,  the  lamer  sca^e  surface  rouahness  usually 
characterized  from  reflectance  measurements  in  terms  of  a  low-wavevector 
distribution.  It  seemed  that  it  would  be  easier  to  observe  splitting  in 
the  surface-plasmon  emission  for  a  given  material  if  the  effects  of  the 
low-wavevector  components  were  minimized.  This  can  be  achieved  in  two 
ways.  For  low-energy  electrons  incident  normal  to  the  surface,  relatively 
few  smal 1 -wavevector  surface  plasmons  are  excited  [1,4].  This  reduces 
the  possible  coupling  between  low-wavevector  surface  plasmons  and  the 
corresponding  low-wavevector  part  of  the  surface-roughness  spectrum,  thus 
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reducing  its  effect  on  the  surface-roughness  coupled  radiative  decay 
spectrum.  In  addition,  the  low-wavevector  components  of  the  surface- 
roughness  spectrum  can  be  made  as  small  as  possible,  which  means  that  on 
a  large  scale  the  film  should  be  as  smooth  as  possible. 

Potassium  is  a  nearly-free-electron  metal  with  a  vol ume-pl asmon 
energy  of  (3.8  ±  0.1)  eV  [5,6].  When  K  films  were  irradiated  with  low- 
energy  electrons  at  near-normal  incidence  [6],  the  surface-pl asmon 
radiation  was  found  to  be  very  sensitive  to  the  temperature  of  the  film 
and  to  the  pressure  in  the  experimental  chamber.  If  the  K  films  were 
not  cooled,  no  surface-pl asmon  radiation  peak  was  observed.  K  has  a 
relatively  low  melting  point  of  62.3°C,  and  it  was  thought  that  the 
smooth  liquid  K  surface  produced  by  electron  beam  heating  could  not 
couple  surface  plasmons  to  the  radiation  field.  The  same  result  was 
obtained  when  the  K  film  was  cooled  to  liquid  nitrogen  temperature,  and 
the  electron  beam  current  was  increased  from  12  to  200  yA. 

Since,  in  order  to  resolve  the  splitting  in  the  surface-plasmon 
radiative  decay  due  to  high-wavevector  components  in  the  surface-roughness 
spectrum,  the  surface  should,  on  a  large  scale,  be  smooth,  it  appeared 
that  K  was  a  good  candidate.  In  addition,  K  has  the  added  advantage  of 
being  free-electron-like  which  provides  a  larger  separation  of  the  volume 
and  surface  pi asmon  energies  than  in  silver.  To  provide  the  high-wave¬ 
vector  components,  this  surface  must  possess  a  periodic  irregularity,  or 
roughness,  with  a  stepping  height  of  the  order  of  a  few  angstroms.  It 
was  felt  that  K  solid  surfaces  allowed  to  solidify  close  to  the  melting 
point  might  be  free  of  stress-induced  large-scale  irregularities  and, 
hence,  show  the  desired  characteristics. 
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Details  of  the  experimental  arrangement  have  been  described  previously 

[6].  A  McPherson  Model  218  scanning  monochromator  was  used  to  measure 

the  spectral  distribution  of  p-polarized  radiation  emitted  from  vacuum- 

evaporated  thick  K  films  bombarbed  by  low-energy  electrons.  A  thick  K 

film  was  vacuum  evaporated  in  situ  directly  onto  a  clean,  shiny  Cu 

substrate  which  was  cooled  by  liquid  nitrogen.  The  evaporator  consisted 

of  a  copper  tube  which  was  attached  to  the  experimental  chamber  through 

a  copper  flange  gasket  and  sealed  off  at  its  other  end.  A  glass  ampoule 

containing  K  was  mounted  in  the  copper  tube.  With  the  experimental 
-8 

system  at  7x10  torr  the  glass  ampoule  was  crushed  by  squeezing  the 

copper  tube.  The  K  was  first  outgassed  for  about  10  minutes  and  then 

evaporated  by  a  heating  coil  wound  around  the  copper  tube.  The  pressure 

rose  as  high  as  10  ^  torr  during  the  evaporation  but  afterwards  dropped 
-8 

rapidly  to  2x10  torr.  1.5-keV  electrons  with  a  beam  current  in  the  10 
-►  20  pA  range  were  incident  on  the  K  film  at  35°  from  the  film  normal, 
and  the  photons  emitted  along  the  normal  were  observed  over  the  wave- 

o 

length  range  from  3000  to  5500  A.  This  range  of  beam  currents  minimized 
heating  of  the  film,  while  still  producing  emission  with  a  reasonable 
intensity.  Following  the  evaporation,  during  which  the  K  film  was 
cooled  with  liquid  nitrogen,  the  coolant  was  removed  and  the  film 
annealed  at  room  temperature.  This  process  of  cooling  the  sample  by 
liquid  nitrogen  and  then  removing  the  coolant  was  repeated.  At  each 
stage  of  this  temperature  cycle,  emission  spectra  were  recorded.  During 
the  second  period  of  annealing,  evidence  was  seen  of  splitting  in  the 
surface-plasmon  emission.  Additional  K  was  then  evaporated  onto  the 
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sample  which  was  now  water  cooled  and  the  emission  spectrum  recorded  as 
a  function  of  time  after  deposition.  Figure  1  shows  the  original  recorder 
trace  obtained  30  minutes  after  the  second  K  deposition.  Figure  2  shows 
the  spectrum  after  correcting  the  spectrum  ir  Fig.  1  for  grating  eff' 

O 

ciency  and  photomul tipi ier  response.  The  peak  at  3250  A  (3.8?  eV)  it 
the  volume-plasmori,  or  transition  radiation,  peak.  When  this  is  sub¬ 
tracted,  peaks  are  seen  at  3980  ft  (3.12  eV)  and  ^ 790  ft  (2.59  eV)  as 
indicated.  These  results  are  representative  of  the  eight  spectra  recorded, 
four  before  and  four  after  the  second  K  deposition.  The  relative  magnitudes 
of  the  peaks  varied  slightly  from  one  spectrum  to  the  next,  the  3980  ft 
peak  generally  growing  with  time  relative  to  the  other  two  peaks.  The 
peak  wavelengths  were  in  the  ranges  (3250  ±  20',  ft,  (3980  t  40)  ft,  and 
(4790  ±  40)  ft  with  no  systematic  variation  with  time  within  these  ranges. 

It  is  thought  that  the  peaks  at  3980  ft  and  47 .!U  ?  are  due  to  surface- 
plasmon  radiative  decay,  on  either  side  of  the  expected  position  of  the 
asymptotic  surface-plasmon  radiation  at  4320  ft. 

The  splitting  of  the  dispersion  relation  for  surface  plasmons  on  a 
rough  solid-vacuum  boundary  for  high  wave' ectors  can  be  expressed  in 
terms  of  the  two  values  of  e, ,  the  real  part  of  the  dielectric  function 
c(cj),  at  which  surface  plasmons  are  generated.  It  has  been  shown  [2] 
that  peaks  occur  at  values  of  given  by  a  and  8  where  a  =  1/B.  For 
reaoy  reference,  c-i  as  a  function  of  wavelength,  taken  from  Raf.  5,  is 
shown  in  Fig.  2.  It  is  seen  that  at  the  3980  ft  peak,  e-j  =  -0.665  =  cr, 
while  at  the  4790  ft  peak,  r.,  =  -1.50  =  6.  Hence,  the  predicted  relation¬ 
ship  between  the  values  of  at  the  peaks,  given  by  a  =  1/8,  is  satisfied 
for  the  peaks  observed  at  3980  ft  and  4790  ft  in  potassium. 


It  has  also  been  shown  [2]  that 
|Ae.j  |  =  | a-S I  <  4a 

where  a  is  a  function  of  the  surface-plasmon  wavevector  and  of  the  high- 

wavevector  components  of  the  surface-roughness  spectrum.  If  a  =  0,  then 

o  =  6  =  -1  and  there  is  no  splitting  of  the  dispersion  relation. 

Figure  2  in  Ref.  2  gives  the  values  of  c-j  at  which  surface-plasmon 

2 

excitation  is  possible.  The  curve  for  a  =0  shows  a  single  peak  at  the 
asymptotic  surface-plasmon  wavelength.  However,  surface  plasmons  gen¬ 
erated  at  a  smooth  surface  cannot  radiate.  The  curves  in  Fig.  2,  Ref. 

2,  for  a  f  0  illustrate  the  expected  splitting  for  surfaces  having  a 
high-wavevector  component  in  the  surface-roughness  spectrum.  For  our 
sample 

| Ae-j  |  =  (-0.67  -  (-1.50)1  =  0-83. 

Hence,  a  ^  0.2. 

A  calculation  of  the  predicted  shape  of  the  surface-plasmon  excita- 
2 

tion  spectrum  for  a  =  0.04  was  performed  as  in  Ref.  2  using  the  published 
values  of  e,  and  e?  f°r  K  given  by  Palmer  and  Schnatterly  [5].  This  is 
shown  in  Fig.  2  for  comparison  with  our  experimental  spectrum. 

Palmer  and  Schnatterly  [5,  Fig.  8]  observed  some  structure  in  an 
absorption  peak  attributed  to  surface-plasmon  excitation  in  K.  This 
structure  was  small  and  was  located  at  2.7  eV  and  2.9  eV,  on  either  side 
of  the  ?.8-eV  surface-plasmon  absorption.  It  is  suggested  that  this 
could  be  splitting  due  to  coupling  to  a  high-wavevector  surface  rough¬ 
ness,  superimposed  on  a  surface-plasmon  peak  due  to  coupling  to  the  low- 
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wavevector  surface  roughness.  From  Fig.  2,  at  2.7  eV  c-j  =  -1.28  =  a, 
while  at  2.9  eV  =  -0.96  =  8-  In  this  case  u  is  only  approximately 
equal  to  1/6-  However,  considering  the  uncertainties  in  the  measured 
quantities,  their  data  are  not  inconsistent  with  splitting  due  to  high 
wavevector  surface  roughness. 
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Fig.  1.  Recorder  trace  of  emission  spectrum  from  K  film  irradiated 
by  1.5-keV  electrons. 

Fig.  2.  Experimental  emission  spectrum  from  K  held  at  near  room 
temperature  and  irradiated  by  1.5-keV  electrons  at  35°  from  the  film 

nornl ;  -  total  measured  intensity,  estimated  transition 

radiation  contribution,  and  -  surface  plasmon  contribution.  Also 

shown  are  the  experimental  c-j  values  from  Ref.  5  (• - •)  and  the 

theoretical  surface-plasmon  spectrum  calculated  as  in  Ref.  2. 
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CHAPTER  III 


Optic al  Properties  of  Polyethylene 

L.  R.  Painter, +  E.  T.  Arakawa,  and  M.  W.  Williams 
Health  and  Safety  Research,  ORNL,  TN  37830 

Abstract 

The  optical  properties  of  polyethylene  in  the  form  of  thin  films 
have  been  determined  for  photon  energies  from  0.5  to  76  eV  from  trans¬ 
mission  measurements.  The  results  for  the  extinction  coefficient,  k, 
have  been  combined  with  calculated  values  in  the  soft  and  hard  x-ray 
regions  and  their  internal  consistency,  and  their  consistency  with  data 
obtained  previously  for  polystyrene,  verified  through  a  sum  rule  calcu¬ 
lation.  It  is  hoped  to  use  these  data  to  calculate  mean  free  paths, 
stopping  powers,  and  "continuous-slowing-down-approximation"  ranges 
for  low-energy  electrons  (MO  eV  to  10  keV)  in  polyethylene. 

Introduction 

Polyethylene  is  used  extensively  for  high-voltage  cable  insulation. 
It  is  also  employed  in  communications  systems,  many  of  which  ar> 
miniaturized  for  use  in  the  space  and  defense  programs.  Hence  an 
understanding  of  the  interactions  of  photons  and  of  high-energy  charged 
particles  with  thin  films  of  polyethylene  is  of  interest.  A  convenient 
way  to  obtain  this  information  is  by  measurements  of  the  optical 
properties  of  thin  film  samples  over  an  energy  range  sufficiently 
large  to  include  all  of  the  oscillator  strength.  From  a  chemical 
viewpoint,  polyethylene  is  the  simplest  synthetic  hydrocarbon  polymer, 
with  a  monomeric  repeat  unit  of -fCH^-CH^]-,  It  has  only  single  bonds 
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and  hence  all  its  valence  electrons  are  designat'd  as  o  electrons. 

We  have  previously  studied  polystyrene^— one  of  the  simplest  of  the 
hydrocarbon  polymers  which  includes  an  aromatic  component  in  the  mono¬ 
meric  unit  and  which  thus  has  both  it  and  o  valence  electrons.  Com¬ 
parison  of  the  optical  properties  of  polyethylene  with  those  for 
polystyrene  extends  our  knowledge  and  understanding  of  the  way  in 
which  electrons  in  hydrocarbon  solids  react  to  electromagnetic  radiation. 

Previous  measurements  of  the  optical  properties  of  polyethylene 
have  been  mainly  restrictea  to  the  IR  and  visible  regions.  In  this 
study,  we  report  the  optical  properties  of  polyethylene  in  the  form 
of  thin  solid  films  over  the  range  of  photon  energies  from  0.5  to 
76  eV.  The  measured  results  for  the  extinction  coefficient,  k,  are  combined 
with  values  calculated  from  those  for  polystyrene, ^ '  to  obtain  a 
combined  k  spectrum  which  covers  substantially  all  of  the  total 
oscillator  strength,  including  the  carbon  core  excitations  starting 
at  eV.  A  Kramers -Kronig  analysis  then  yields  the  refractive 
index,  n,  from  0.5  to  76  eV,  and  a  sum-rule  calculation  verifies  the 
internal  consistency  of  these  optical  properties.  Comparisons  ace 
drawn  between  the  measured  optical  properties  of  polyethylene  and 
polystyrene.  It  is  hoped  that  differential  inverse  mean  free  paths 
(DIMFP's)  can  be  calculated  for  the  excitation  of  inner-shell  electrons 
in  polyethylene.  These  combined  with  DIMFP’s  from  the  model  insulator 
theory  would  be  used  to  calculate  electron  mean  free  paths  and  stopping 
powers  for  comparison  with  the  Bethe-Bloch  theory  at  high  energies.  If 
these  comparisons  are  satisfactory  a  complete  set  of  tabulations  for 
polyethylene,  as  has  been  done  previously  for  polystyrene,^'^  will 
be  prepared. 
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Experiment a  J  Techniques 

(4) 

Transmission  techniques  were  used  to  obtain  values  of  the  extinction 

coefficient,  k,  of  polyethylene  films  over  essentially  the  whole  energy 

range  from  0.5  to  76  eV.  Up  to  6.6  eV,  transmission  measurements  were 

(5) 

made  on  free-standing,  0.5  mil.  Glad  Wrap  '  using  a  Cary  spectrometer 
and  a  Seya-Namioka  monochromator.  In  this  case,  the  transmittance,  T, 
yielded  absolute  values  of  k  since  the  thickness,  t,  of  the  film  was 
k nown,  and 

T  9  -4nkt 

T  =  T~  =  (1  -  R;  e  1  . 

o 

I  was  the  intensity  of  the  radiation  incident  normally  on  the  poly¬ 
ethylene,  while  I  was  the  intensity  transmitted  by  the  film.  R  was 
the  normal  incidence  reflectance  for  the  polyethylene-air  interface 
and  was  calculated  from  the  n  values  measured  on  evaporated  polyethylene 
films,  using  methods  described  later.  Below  about  4.9  eV  (k  =  10"^) 
k  was  too  small  to  be  measured  using  a  0.5  mil  thick  sample.  The  use 
of  thicker  free-standing  polyethylene  samples  (0.9  mil  Glad  Wrap  and 
4  mil  tialgene  polyethylene)  showed  weak  absorption  below  n/|  eV  but  no 
measurable  absorption  from  n/|  eV  to  n4.5  eV.  Data  above  6.6  eV  were 
obtained  on  evaporated  films  of  polyethylene  prepared  by  vacuum  deposi¬ 
tion  from  Dow  Corning^  polyethylene  pellets.  However,  since  it  proved 
difficult  to  obtain  vacuum-evaporated,  free-standing  polyethylene  films 
suitable  for  transmission  measurements,  techniques  were  used  which 
allowed  the  polyethylene  to  be  supported  by  a  substrate.  For  use  in 
the  7.7  to  10.6  eV  energy  range,  polyethylene  was  evaporated  over  half 
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of  a  MgF^  disc.  Transmission  measurements  were  made  by  recording  I 
through  the  MgF£  disc  alone  and  I  through  the  KgF^  plus  polyethylene 
film,  using  a  Seya-Namioka  monochromator.  Two  separate  films  were 
used  with  thicknesses  of  460  and  504  A  as  determined  by  interferometric 
and  by  quartz  crystal  measurements.  These  measurements,  limited  by 
the  absorption  edge  of  the  MgF^  substrate  at  MO. 6  eV,  also  allowed  an 
absolute  determination  to  be  made  of  the  k  values.  Above  10.6  eV  the 
fluorescent  substrate  method^  was  used.  A  glass  slide  was  coated  with 
a  thick  layer  of  p-terphenyl ,  which  is  fluorescent.  Polyethylene  was 
then  vacuum  evaporated  over  half  of  the  p-t.erpheny'' .  Radiation  was 
incident  normally  on  this  sample  and  the  fluorescent  radiation  from 
the  p-terphenyl  monitored.  For  the  uncovered  p-terphenyl  the  intensity 
of  the  fluorescent  radiation  was  «I  ,  while  for  the  polyethylene-covered 
p-terphenyl  it  was  <*1.  Since  p-terphenyl  fluoresces  in  the  visible,  the 
intensity  of  the  fluorescent  radiation  was  not  attenuated  on  passing  out 
through  the  glass  slide  and  hence  the  measured  intensities  from  the  two 
halves  of  the  sample  were  subject  to  the  same  proportionality  constant 
and  their  ratio  gave  the  transmittance,  T,  of  the  polyethylene.  From 
10.6  to  76  eV  observations  were  made  using  a  grazing-incidence  mono¬ 
chromator.  From  10.6  to  M35.5  eV  where  k  is  relatively  large,  a  thin 
polyethylene  film  was  used.  Its  thickness  was  estimated  as  M00  A  by 
normalizing  k  to  the  experimentally  determined  absolute  value  at  10.6  eV. 
This  process  was  repeated  from  35.5  to  76  eV,  using  a  thicker  film  which 
was  estimated  to  be  Mi90  A  thick,  by  normalizing  k  at  35.5  eV.  The 
values  of  k  obtained  by  the  fluorescent  substrate  method  that  are 
presented  in  this  paper  were  confirmed  on  a  number  of  other  films  with 
similar  thicknesses  to  those  described  above. 
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Values  of  the  real  part,  n,  of  the  complex  refractive  index, 
n  =  n  +  ik,  were  also  determined  experimentally  in  the  energy  range 
from  1.77  to  8.4  eV.  From  1.77  to  5.8  eV,  n  was  determined  from  critical 
angle  measurements. Polyethylene  was  vacuum  evaporated  onto  the 
plane  face  of  a  sapphire  semicyl  ir.der .  Radiation  from  a  beya-Mamioka 
monochromator  was  i nc i den t  along  a  radius  of  the  semicylinder,  reflected 
from  the  sapphire-polyethylene  interface,  and  exited  along  a  radius. 

Scans  of  reflectance  v_s  angle  of  incidence  at  the  sapphire-polyethylene 
interface  were  recorded.  For  photon  energies  where  the  refractive  index, 
r  ,  of  the  material  of  the  semicylinder  is  greater  than  that  of  the  poly¬ 
ethylene,  n ,  sin<c  =  n/r,  where  the  critical  angle,  ec,  is  defined  as  the 
angle  which  marks  the  onset  of  total  internal  reflection.  Thus  if  ri  is 
known  and  0C  measured,  n  can  be  calculated.  In  this  case  values  of  n$ 
were  obtained  before  the  polyethylene  was  deposited  on  the  semicylinder 
from  scans  of  reflectance  v£  angle  of  incidence  at  the  plane  sapphire- 
air  interface.  From  7.3  to  8.4  eV  reflectance  vs  angle  of  incidence 
scans  were  ootained  at  the  CaF  -polyethylene  interface  for  polyethylene 
on  a  CaF;  semicylinder.  In  this  case  critical  angles  were  not  observed, 
and  the  reflectance  vs  0  data  were  least  squares  fit  to  Fresnel’s  equa¬ 
tions.  Interference  fringes  were  also  analyzed  for  the  polyethylene  films 
on  the  MgF^  discs  to  yield  both  n  and  the  film  thicknesses.  The  values 
of  n  obtained  from  each  of  these  techniques  agree  well  with  the  two 
values  available  in  the  1 iterature. 
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Results  and  P 

The  experimental  results  obtained  f  .  k  •  •  .(  a^d  76  eV 

are  presented  in  Fig.  1  on  a  logarithmic  pi.  The  best  smooth  curve 

through  the  experimental  points  is  also  shown  in  Fig.  1.  The  same  data  is 

presented  in  Fig.  2  on  a  linear  plot  from  0.5  to  76  eV.  Also  presented 

in  Fig.  2  are  the  n  values  determined  experimentally  in  the  present  study 

and  from  the  1 iterature. A  Kramers-Kronig  analysis  was  performed  on 

these  data  in  order  to  obtain  values  of  n  for  polyethylene  over  the  whole 

experimental  range  from  0.5  to  76  eV.  The  k  values  from  o,30  eV  to  76  eV 

were  extrapolated  as  shown  in  Fig.  1.  Since  absorption  measurements  above 

76  eV  were  unavailable  for  polyethylene,  the  contribution  due  to  the 

carbon  Is  electrons  was  estimated  from  data  for  polystyrene.^  It  was 

assumed  that  in  hydrocarbons,  which  contain  only  carbon  and  hydrogen,  the 

experimental  k  associated  with  the  C  Is  electrons  is  proportional  to  the 

number  of  C  Is  electrons  per  unit  volume.  The  contribution  to  k  for 

polyethylene  due  to  core  electrons  was  calculated  on  this  assumption 

above  the  carbon  K-shell  absorption  edge  from  experimental  data  for  poly- 

-3  -3 

styrene,  assuming  densities  of  1.05  g.cm  for  polystyrene  and  0.92  g.cm 
for  polyethylene.  The  calculated  core  electron  contribution  to  k  was 
then  added  to  the  extrapolated  valence  electron  contribution  as  shown 
in  Fig.  1  to  give  the  total  estimated  k-spectrum  for  polyethylene. 

The  refractive  index,  n,  as  a  function  of  photon  energy,  E,  was 
determined  from  the  k-spectrum  by  means  of  the  Kramers-Kronig  relation 


n(E)  -  1 


2  f"  E'ME’) 
o  E^TT2 


dE  ’  . 
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In  the  region  below  V!  eV,  where  photoabsorption  due  to  molecular 
vibrations  occurs,  the  k  values  were  so  small  that  the  contributions 
to  the  integrals  of  E’k(E'),  in  this  low-energy  region,  were  neglected. 

At  high  energies  the  cut-off  for  the  integral  was  chosen  so  that  the 
values  of  n  obtained  from  the  Kramers-Kronig  analysis  agreed  with  tiiose 
determined  e.>pe  imentally  in  the  1.77  to  3.4  eV  range.  The  values  of  n 
obtained  for  polyethylene  by  this  analysis,  from  0.5  to  76  eV ,  are  shown 
in  Tig.  2.  Then  assuming  n  =  n  +  ik  and  =  e  =  +  ic^.  The  n  and  k 

values  were  used  to  calculate  the  real  and  imaginary  parts  f-|  and  ?  of 
the  complex  dielectric  function,  e,  and  the  energy  loss  function  In(-l/e). 
These  quantities  are  shown  in  Figs.  3  and  4,  respectively.  The  estimated 
error  on  the  experimental  k  values  is  6. +6%  and  on  the  experimental  n  values 
is  £±0.5%.  The  possible  errors  on  the  other  quantities  depend  on  the 
validity  of  the  Kramers-Kronig  analysis. 

A  check  was  performed  on  the  internal  consistency  of  the  optical 
properties  derived  for  polyethylene  and  on  their  consistency  with  data 
obtained  previously  for  polystyrene.  For  e2  and  Im(-l/c'),  there  exist 
the  well-known  f  sum  rules^^ 

E 

ME)  =  |  -  C  Fie2(E')JE,->  N,  E  -  -  , 

47rn„e  dT  J 
o  o 

and 

N2<e’  '  ’  dE'  *N-E*'”' 
ft  n  ' 
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where  n_  is  the  molecular  density  and  N  is  the  total  number  of  electrons 
in  a  molecule.  Taking  a  monomeric  unit  — of  polyethylene  as  a 
molecular  unit,  N  is  16.  The  calculated  values  of  N-j(E)  and  ^(E)  as 
functions  of  photon  energy,  E,  are  shown  in  Fig.  5.  Both  N^(E)  and 
^(E)  are  seen  to  saturate  at  16  electrons  per  molecular  unit,  indicating 
interna1  consistency  of  the  optical  constants  determined  for  polyethylene. 
Since  the  core  contribution  to  k  was  calculated  from  polystyrene  data, 
this  also  shows  the  validity  of  our  assumptions  concerning  the  general 
behavior  of  core  electrons  in  hydrocarbons  containing  only  carbon  and 
hydrogen . 

Further  evaluations  of  N( E)  were  made  on  the  k  values  due  to  the 
valence  and  core  excitations,  separately.  The  results  for  NV(E)  and 
NC(E),  the  contributions  to  N-j(E)  or  N^(E)  above  the  carbon  K-shell 
absorption  edge  from  the  valence  and  core  electron  excitations,  respec¬ 
tively,  are  also  shown  in  Fig.  5.  A  redistribution  of  electron  numbers 
corresponding  to  0.5  electrons  is  found  between  the  valence  and  core 
excitations.  NV(E)  and  NC(E)  should  tend  to  12  and  4,  the  numbers  of 
valence  and  core  electrons,  respectively,  in  the  limit  E  -*■  °°,  if  there 
were  no  oscillator  strength  coupling  between  them.  A  similar  result  was 
found  previously  for  polystyrene.^  As  for  polystyrene,  it  should  also 
be  noted  that  the  effective  numbers  of  electrons  as  plotted  in  Fig.  5 
as  functions  of  photon  energy  E  do  not  saturate  until  a  few  thousand  eV. 
This  is  so  even  for  the  valence  electrons  shewing  that  even  at  these 
high  energies  the  contribution  of  valence  excitations  to  the  oscillator 
strength  is  not  negligible.  For  both  polyethylene  and  polystyrene, 
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these  results  demonstrate  that  sum  rule  calculations  can  only  be  applied 
to  the  total  number  of  electrons  in  the  system.  Partial  sum  rules, 
employed  in  some  other  situations,  cannot  be  used  for  hydrocarbons. 

Comparison  of  the  plots  of  n,  c,  and  lm(-l/e)  for  polyethylene  with 
the  corresponding  variations  with  energy  for  polystyrene^  shows  striking 
similarities.  The  only  fundamental  difference  is  the  sharp  structure 
seen  in  polystyrene  below  n/|0  eV  associated  with  the  tt  electrons  which 
is  not  present  in  polyethylene.  Thus  the  general  conclusions^  arrived 
at  from  analysis  of  sum  rules  for  polystyrene  concerning  the  interaction 
of  electromagnetic  radiation  with  the  electrons  in  a  hydrocarbon  solid 
over  essentially  the  whole  energy  range,  can  be  applied  to  polyethylene. 
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Figure  Captions 

Fig.  1.  Extinction  coefficient,  k,  for  polyethylene  as  a  function 
of  photon  energy.  The  calculated  values  were  obtained  from 
experimental  values  for  polystyrene.^ 

Fig.  2.  Optical  properties  of  polyethylene  over  the  energy  range 
of  the  experimental  measurements.  The  smooth  curve  for  n  was 
derived  from  a  Kramers-Kronig  analysis  of  the  k  values.  The 
experimental  values  of  n  are:  ORNL,  and  from  Ref.  IO. 

Fig.  3.  Dielectric  functions  of  polyethylene  as  a  function  of  photon 
energy. 

Fig.  4.  Energy  loss  function  of  polyethylene  as  a  function  of  photon 
energy. 

Fig.  5.  Effective  number  of  electrons  N^(E)  and  ( E )  per  monomeric 

unit  of  polyethylene  obtained  from  the  imaginary  part  of  the 

A* 

dielectric  function  e  =  +  ie^,  and  the  energy-loss  function 

lm(-l/c),  respectively,  as  functions  of  photon  energy  E.  NV(E) 
and  N  (E)  represent  the  contributions  from  the  valence  and  core 
electron  excitations,  respectively. 
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CHAPTER  IV 


Inelastic  Interactions  of  Electrons  with  Polystyrene:  Calculations 
of  Mean  Free  Paths,  Stopping  Powers,  and  CSDA  Ranges 

A  theoretical  description  of  the  inelastic  interactions  of  electrons 
with  solid  polystyrene  has  been  presented.  The  response  of  the  valence 
electrons  to  energy  and  momentum  transfers  was  determined  by  a  model 
insulator  theory;  carbon  K-shell  ionization  cross  sections  were  derived 
from  atomic,  generalized  oscillator  strengths.  Contributions  to  the 
inverse  mean  free  paths  and  stopping  powers  due  to  these  two  excitation 
processes  were  derived  and  tabulated  for  incident  electrons  with  energies 
from  10  eV  to  10  keV.  Electron  ranges  in  the  continuous-slowing-down 
approximation  were  calculated  and  tabulated  for  electrons  with  energies 
from  15  eV  to  10  keV.  Complete  descriptions  of  this  calculational  method 
and  these  results  may  be  found  in  the  following  references: 

1.  J.  C.  Ashley,  C.  J.  Tung,  R.  H.  Ritchie,  and  V.  E.  Anderson, 
RADC-TR-78-32  (February,  1978),  A057881 . 

2.  J.  C.  Ashley,  C.  J.  Tung,  and  R.  H.  Ritchie,  IEEE  Trans.  Nucl. 
Sci.  NS-25,  1566-1570  (December,  1978). 
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CHAPTER  V-A 


The  ATR  Method  with  Focused  Light--Appl ication  to  Guided  Waves 

on  a  Crating 

The  attenuated  total  reflection  (ATR)  method  has  been  modified  by 
using  focused  light.  The  reflected  light  cone  contains  information 
about  a  large  portion  of  wavevector  space  instead  of  only  one  point 
as  in  the  conventional  ATR  method.  The  dispersion  relation  of  guided 
modes  can  be  found  directly  from  the  image  of  the  light  cone  on  a 
screen.  When  applied  to  guided  waves  on  a  gratinq,  the  dispersion 
relation  of  the  waves  with  "band  gaps"  and  discontinuities  can  be 
seen  on  the  screen.  This  application,  with  further  discussion  of 
the  ATR  method  using  focused  light,  has  been  published.  See: 

E.  Kretschmann,  Opt.  Commun.  26,  41-44  (July,  1°78). 
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CHAPTER  V-B 


Refractive  Index  of  LiF  Films  as  a  Function  of  Time* 


T.  L.  Haltom,  E.  T.  Arakawa,  M.  W.  Williams  and  E.  Kretschmann 
Oak  Ridge  National  Laboratory,  Health  and  Safety  Research  Division, 
Oak  Ridge,  Tennessee  37830,  USA 


★ 

Research  sponsored  jointly  by  the  Southern  College  University 
Union  the  Deputy  for  Electronic  Technology,  Air  Force  Systems  Command, 
under  Interagency  Agreement  DOE  No.  40-226-70  and  the  Division  of 
Biomedical  and  Environmental  Research,  U.  S.  Department  of  Energy,  under 
contract  W-7405-eng-26  with  the  Union  Carbide  Corporation. 

Also  Millsaps  College,  Jackson,  MS  39202. 

Present  address:  Schaeferstr.  28,  2  Hamburg  6,  West  Germany 


117 


Abstract 


The  optical  properties  and  thickness  of  evaporated  Lit7  films  on 
silver  or  aluminum  were  determined  by  attenuated  total  reflection  tech¬ 
niques  as  a  function  of  time  after  deposition.  With  the  L i F  films 
maintained  in  vacuum,  the  refractive  index  increased  from  1.3  to  "vl.4, 
or  close  to  the  crystalline  value,  in  the  course  of  a  few  days.  This 
implies  that  the  LiF,  loosely  packed  when  first  deposited,  anneals  over 
a  period  of  time  to  the  more  closely  packed  crystalline  configuration 
Tne  LiF  films  were  found  to  be  ootically  isotropic  in  contras'*  to  earlier 
reports  of  anisotropy. 

Int rod  uction 

Attenuated  total  reflection  (ATR)  techniques  employing  various 
experimental  configurations  have  been  reviewed  by  Otto. '  viore  specifi¬ 
cally,  such  techniques  have  been  used  to  study  light  modes  in  evaporated 

films  of  LiF  and  the  results  analyzed  in  terms  of  the  optical  properties 

2  **4 

and  thickness  of  the  films.  However,  the  behavior  of  these  properties 
with  time  after  deposition  has  net  previously  been  studied  for  LiF 
films.  The  changes  in  the  refractive  index,  n,  reported  here  foe  LiF 
indicate  that  changes  take  place  in  the  film  structure  after  deposition. 
These  results  become  significant  when,  for  example,  an  optical  system 
involving  reflection  or  antireflection  coatings  must  be  designed  for 
maximum  efficiency.  It  is  then  necessary  to  know  the  correlation  between 
the  properties  of  the  coating  as  deposited  and  in  its  final  state. 
Deterioration  of  coatings  on  exposure  to  certain  anticipated  environ¬ 
mental  conditions  could  also  be  studied  experimentally  by  this  technique. 
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Ex r e r i mental  Ts c h rru ] a c_s 

The  experimental  arrangement  was  similar  to  that  used  previously  by 

A 

Hornauer  a>  d  Raether'*  tor  LiT.  A  transparent  semicyl inder  was  mounted 
in  ,  ref  1  ectometer  in  ^ucr  a  /ay  that  a  narrow  bean  of  monochromatic 
light  entered  the  semi cy ) i nder  along  a  radius  and  hence  normal  to  the 
curved  surface,  was  reflected  at  the  plane  semicyl inder-vacuum  interface, 
and  exited  from  the  semi cvl i nder  along  a  radius.  The  incident  liqht  was 
P7ane  polarized  /sing  a  Polaroid  filter,  and  the  reflected  light  was 

-6 

detected  with  a  photomultiplier.  The  pressure  in  the  system  was  10 
torr. 

A  thi  metallic  layer  wa :  vacuum  evaporated  in  sjtu  onto  the  plane 
'■‘ace  of  the  semicylinder.  A  scan  of  reflectance  vs  angle  of  incidence 
for  p-pol irized  light  showed  a  strong  dip  in  reflectance  at  an  angle  of 
incidence  just  above  the  critical  angle  due  to  the  generation  of  surface 
plasmons  at  the  metal -vacuum  interface.  A  thick  Li F  film  was  then 
deposited  by  /acuum  evaporation  over  the  thin  metal  film.  Scans  of 
r  flectance  vs  anqle  of  incidence,  recorded  indi vidual ly  for  p-  and  s- 
polarizei  light,  then  showed  several  strong  dips  at  angles  above  the 
critical  angle.  There  are  two  types  of  resonances  possible  in  the 
semicyl inder/metal /Li F  system:  surface  plasmons  and  ontical  guided 
waves.  The  surface  plasmon  mode,  which  now  occurs  at  the  metal -Li F 
interface  is  shifted  to  a  larger  angle  of  incidence  than  its  value  for 
the  metal-vacuum  interface.  In  both  systems  the  dip  corresponding  to 
the  surface  plasmon  mode  shows  broadening  resulting  mainly  from  the 
energy  absorption  associated  with  the  non -zero  value  o'  the  extinction 
coefficient  of  the  metal.  Several  sharper  dips  in  reflectance  showinq 
slight  instrumental  broadening  only  and  corresponding  to  resonant 


119 


conditions  for  optical  guided  modes  of  different  orders,  occur  between 
the  critical  ang1^  and  the  surface  plasmon  resonance.  The  thicker  the 
LiF  film,  the  greater  the  number  of  guided  modes  that  can  be  supported. 

With  the  system  maintained  under  the  'vl 0 torr  vacuum,  reflectance  vs 
angle  of  incidence  scans  were  repeated  at  different  times  after  LiF  film 
deposition.  It  was  found  that  the  characteristic  reflectance  minima 
always  shifted  to  larger  angles  of  incidence  with  increasing  age  of  the 
LiF  film,  each  minimum  going  asymptotically  with  time  to  its  final 
angular  position.  Observations  were  repeated  until  the  characteri Stic 
minima  were  effectively  no  longer  changing  position  with  time.  Repre¬ 
sentative  reflectance  vs  angle  of  incidence  scans  obtained  using  p- 
polarized  light  incident  on  a  semicyl inder/metal /Li F  system  are  shown  in 
Fig.  1.  In  this  case  the  semicylinder  was  of  fused  silica,  the  metal 
was  Ag,  the  LiF  was  calculated  to  have  a  thickness  of  855  nm,  and  the 
wavelength  of  the  light  was  500  nm.  The  scan  for  time  zero  was  recorded 
immediately  after  deposition  of  the  LiF  and  shows  both  optical  guided 
modes  and  the  surface  plasmon  mode  (labelled  m  =  0).  For  the  semicylinder/ 
Ag  film  alone  the  surface  plasmon  mode  was  excited  at  an  angle  of  incidence 
of  ^46",  whereas  it  appeared  initially  for  this  semicyl inder/Aq/Li F 
system  at  V75°.  The  scans  recorded  after  two  and  after  thirteen  days 
illustrate  that  for  a  particular  optical  mode  (e.g.,  labelled  m  =  1) 
the  shift  to  higher  values  of  the  observed  angle  of  incidence  for  resonance 
is  most  rapid  in  the  first  few  days. 

Observations  were  obtained  on  several  systems,  for  various  periods 
of  time  lasting  from  one  to  26  days.  In  each  case  a  Suprasil  (high- 
purity  fused  silica)  semicylinder  was  used,  and  the  metal  layer  was  Aa 


or  A1 .  The  thin  metal  films  ranged  in  thickness  from  20  to  50  mn  for  Ag 
and  from  15  to  30  nm  for  A1 .  The  L i F  films  ranged  from  855  to  2,100  nm 
thick. 

Analysis 

The  reflectance  vs  angle  of  incidence  scans  obtained  on  the  thin 
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metal  films  were  analyzed  by  a  method  used  previously  by  Kretschmann. 
Through  the  use  of  approximate  relationships  involving  the  angular 
position  of  the  surface  plasmon  reflectance  minimum,  its  depth,  and  its 
half-width  the  real  and  imaginary  parts  of  the  complex  dielectric  constant 
and  the  thickness  of  each  metal  film  were  obtained.  These  quantities 
are  required  in  a  proper  interpretation  of  the  bulk  properties  of  the 
L i F  from  the  positions  of  the  guided  mode  resonances. 

The  guided  mode  resonances  in  the  Lif  were  analyzed  using  relation- 
ships,  given  by  Hornauer  and  Raether,  which  allow  for  the  possibility 
of  the  LiF  being  optically  anisotropic.  Small  corrections,  6@0,  were 
first  madeb  to  the  observed  angles  of  incidence  at  which  resonance 
occurred,  p,  in  order  to  obtain  the  correct  values,  6  ,  in  the  absence 
of  couplino  between  the  L i F  and  the  semicylinder.  It  is  the  thickness 
of  the  metal  film  which  determines  the  extent  of  this  coupling.  We 
chose  the  metal  thickness  to  produce  good  coupling  with  only  slight 
distortion  (60Q  <<  6  )  of  the  guided  modes  in  the  LiF.  In  Hornauer  and 
Raether' s  experiment  this  correction  was  not  necessary  as  their  Ag  films 
were  '^80  nm  thick  and  the  LiF  was  almost  completely  decoupled  from  their 
prism.  A  least  squares  fitting  procedure  was  used  to  fit  simultaneously 
the  corrected  angular  positions  of  all  of  the  observed  p-  and  s-guided 
mode  resonances,  for  a  given  film  and  wavelength,  to  Eqs.  (3)  in  Ref. 
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(4)  to  yield  the  film  thickness  D,  and  nj_  and  nj , ,  the  refractive 
indices  of  the  LiF  for  light  polarized  perpendicular  and  parallel, 
respectively,  to  the  plane  of  incidence.  The  only  mode  not  included  in 
this  analysis  was  the  surface  plasmon  mode,  which  is  associated  with  the 
metal-LiF  interface  rather  than  with  the  bulk  properties  of  the  LiF 
film. 

The  results  of  our  observations  on  several  systems  are  shown  in 
Figs.  2  and  3.  For  all  the  films  studied  the  LiF  was  essentially  opti¬ 
cally  isotropic,  any  deviation  from  isotropicity  being  <0.2%.  Hence  the 
calculated  refractive  index,  n,  of  the  LiF  is  graphed  as  a  function  of 
time,  t,  after  LiF  deposition.  In  Fig.  2  the  values  of  n  vs  t  for  a 
wavelength  of  500  nm  correspond  to  the  reflectance  vs  6  data  shown  in 

Fig.  1.  Also  shown  are  the  n  vs  t  data  for  three  other  wavelengths. 

Figure  3  snows  data,  obtained  at  X  =  500  nm,  for  several  Suprasi 1/Ag/LiF 
and  Suprasi 1 / A1 /Li F  systems.  For  each  sample  the  value  of  n  was  VI.  3 
when  the  film  was  deposited  and  increased  to  VI.4  after  several  days. 

Also  shown  in  Fig.  3  is  the  LiF  thickness,  D  vs^  t,  for  one  of  the  samples. 

In  each  system  a  decrease  in  LiF  thickness  of  V|%  was  seen  in  the  first 
one  to  two  days  after  deposition,  after  which  the  thickness  was  constant 
with  time  to  vtO.2%.  As  a  check  on  the  results  presented,  a  thick  LiF 
film  was  deposited  directly  onto  the  plane  face  of  a  semicylinder  and 
the  refractive  index  of  the  LiF  calculated  from  observations  of  the 
critical  angle  at  the  semicylinder/LiF  interface.  The  variation  of 
refractive  index  with  time  after  deposition  for  this  system  is  also 
shown  in  Fig.  3  and  is  seen  to  be  consistent  with  data  obtained  on  the 
semicyl inder/metal /Li F  systems. 


For  a  single  resonance,  uncertainties  due  to  instrumental  factors, 
including  some  slight  deviation  from  the  described  geometry,  account  for 
an  uncertainty  in  n  of  ±0.002  or  <0.2%.  However,  since  each  n  value 
results  from  a  least  squares  fit  to  the  observed  angular  positions  of 
six  or  more  resonances,  the  estimated  error  in  a  given  n  value  is  <±0.001. 
Discussion 

Any  change  in  the  angular  position  of  a  given  guided  mode  resonance 
implies  a  change  in  the  optical  path  within  the  film  caused  by  a  change 
in  the  refractive  index  and/or  the  thickness  of  the  Li F  film.  We  have 
demonstrated  that  for  LiF  films  maintained  in  vacuum,  changes  in  the 
optical  path  within  the  film  result  mainly  from  a  change  in  the  refrac¬ 
tive  index.  The  refractive  index  increases  from  %1.3  to  VI. 4,  or  close 
to  the  crystalline  value,  in  the  course  of  a  few  days.  This  implies 
that  the  LiF,  loosely  packed  when  first  deposited,  anneals  over  a  period 

of  time  to  the  more  closely  packed  crystalline  configuration.  We  saw  no 

4 

evidence  of  the  V|%  anisotropy  reported  by  Hornauer  and  Raether.  Any 
anisotropy  present  in  our  films  was  <0.2%  as  long  as  the  LiF  was  main¬ 
tained  in  vacuum.  However,  when  observations  were  made  on  a  semicyl inder/ 
Ag/LiF  system  exposed  to  the  atmosphere  and  the  reflectance  vs  0  scans 
analyzed  in  the  same  way  as  the  data  presented  in  this  paper,  the 
refractive  index  of  the  LiF  appeared  to  show  a  small  anisotropy^  similar 
to  that  seen  by  Hornauer  and  Raether.  However,  these  data  could  be 
explained  more  satisfactorily  in  terms  of  a  thin  barrier  layer  with  an 
index  of  refraction  a,  1.15  separating  an  optically  isotropic  LiF  film 
from  the  Ag  substrate.  All  of  the  resonances,  including  the  surface 
plasmon  mode,  are  included  in  this  analysis,  whereas  the  analysis 
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leading  to  the  ^1%  anisotropy  does  not  include  the  surface  plasmon  mode. 
Our  Ag/LiF  samples  exposed  to  the  atmosphere  eventually  peeled,  or 
flaked,  off,  whereas  those  maintained  in  vacuum  remained  stable  in 
appearance.  The  intermediate  layer  apparently  originates  from  the  Ag 

g 

films.  It  has  been  found  that  when  silver  films  were  annealed  in 
oxygen-containing  atmospheres,  surface  agglomeration  occurred  rapidly; 
whereas  when  annealed  in  other  selected  atmospheres,  they  were  stable. 

The  time  required  to  approach  the  crystalline  refractive  index  is 
found  to  vary  among  the  different  samples,  and  the  factor  governing  the 
rate  of  change  of  n  is  not  known.  In  each  case,  however,  the  change  in 
n  with  time  is  most  rapid  immediately  after  film  deposition,  decreasing 
with  increasing  time  as  n  approaches  the  crystalline  value.  Differences 
between  systems  presumably  are  inherent  to  the  samples. 

For  comparison  with  the  values  of  refractive  index  that  we  report 

4 

here,  Hornauer  and  Raether  gave  the  following  results  for  two  different 
films : 


(i) 

=  1.3086, 

n„  =  1.3175, 

D  =  3093  nm 

( i  i )  nx 

=  1.3090, 

n  |  (  =  1 .3160, 

D  =  3143  nm 

3  2 

Otto  and  Sohler  gave  n  =  1.31  for  D  =  3000  nm,  while  Holst  and  Raether 

gave  n  =  1.22  for  D  <  20  nm  and  n  =  1.40  for  D  >  50  nm.  There  is  no 
mention  in  any  of  these  reports  of  any  variation  of  n  with  time. 

Our  results  are  significant  when  an  optical  system  must  be  designed 
for  maximum  efficiency  under  specific  conditions.  Dielectric  layers 
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which  act  as  reflection  or  antireflection  coatings  are  frequently  used 
on  optical  components,  and  the  thickness  desired  for  a  given  function 
depends  on  the  final  value  of  the  refractive  index  of  the  material  of 
the  layer.  Thus  it  is  necessary  to  know  the  correlation  between  the 
properties  as  deposited  and  the  final  properties  of  the  layer.  In  most 
applications  the  exact  time  taken  to  reach  the  final  crystalline-like, 
or  annealed  conditions,  would  be  unimportant.  It  is  sufficient  to  know 
that  an  approximately  1 0%  increase  occurs  in  n,  from  M.3  to  ^1.4,  after 
film  deposition.  Thus  if  films  are  monitored  during  deposition  by 
interferometric  techniques,  subsequent  changes  in  n  should  be  taken  into 
account  when  choosing  the  initial  film  thickness  for  a  particular  desired 
property. 
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Figure  Captions 

F;g.  1.  Measured  reflectance  as  a  function  of  angle  of  incidence, 
on  a  semicyl inder/Ag/LiF  system,  recorded  for  p-polarized  light  of  a 
wavelength  of  500  nm.  The  characteristic  minima  are  seen  to  shift  to 
larger  angles  with  increasing  time  after  Li F  deposition. 

Fig.  2.  The  refractive  index,  n,  of  one  LiF  film  as  unction  of 
time  after  deposition  and  of  the  wavelength,  X,  of  the  incident  light. 

Fig.  3.  The  refractive  index,  n,  for  several  LiF  samples  as  a 
Tunction  of  time  after  deposition,  'or  a  wavelength  of  50C  nm. 

- Suprasi 1 /Ag/Li F  systems.  - Suprasil/Al/LiF  systems. 

....  Suprasil/LiF  system.  The  thickness  of  the  LiF  film  is  recorded 
or.  each  curve.  Also  shown  is  the  thickness  as  a  function  of  time  o^  one 
'.iF  samole. 
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CHAPTER  VI 

A  MODEL  CALCULATION  OF  CROSS  SECTIONS  FOR  INNER-SHELL 
IONIZATION  BY  ELECTRONS 

J.  C.  Ashley  and  R.  H.  Ritchie 
Health  and  Safety  Research  Division 

Inverse  mean  free  paths  (or  cross  sections)  for  inner-shell  ionization 
of  atoms  by  electrons  form  part  of  the  basic  input  for  detailed  electron 
transport  calculations.  Thus  it  is  important  to  have  a  simple,  reliable 
method  for  obtaining  such  cross  sections  for  a  wide  variety  of  materials. 

We  propose  a  scheme  for  predicting  cross  sections  for  inner-shell  ionization 
based  on  binary  encounter  theory  and  experimental  photoelectric  cross 
sections.  Background  material  for  this  calculation  can  be  found  in  our 
earlier  study  entitled  "Classical  Binary  Collision  Cross  Sections  of 
Atomic  Systems"  (see  Ref.  1). 


The  fori' 


g(e)  *  A0  +  A,?.  +  A ztL  +  A3e*  .3; 

is  convenient  and  should  be  flexible  enough  to  fit  ex,,erv)entu  data  ,.er 
a  fairly  .vide  range  of  energy.  In  the  limit  -0,  Eq.  (2)  becomes 


Jjj  _  2° 

o(£  "  3tr 


M) 


Note  that  the  integral  -g—  dL  Mould  diverge  if  g(t)  conta'nec 

a  power  of  e  higher  than  the  third.  If  terms  with  negative  powt-'S 
had  been  included,  the  integral  of  Eq.  (2)  over  t  would  become  very 
cumbersome  to  obtain  in  analytical  form. 

The  constants  required  to  specify  gtc),  an.;  thus  the  optical  oscillator 
strength,  Eq.  (4),  an  be  obtained  from  data  on  photoelectric  cross 
sections.  (£  ,  through  the  relation 


dft> 


_ Op(c) 
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(5) 


where  aQ  is  the  Bohr  radius  and  1  is  the  fine  structure  constant,  or 

[eV]  =  9.  //  *  10' 9 crp  [barns]  .  (£,) 


The  compilations  of  photoabsorption  cross  sections  produced  by  Biggs 
2  3 

and  Lighthill  and  by  Storm  and  Israel  are  useful  for  fitting  purposes. 
Biggs  and  Lighthill  give  analytical  fits  to  experimental  data  for 
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a  wide  range  of  elements  (1  Z  100).  Unfortunately  their  fits  have 
the  form  rr  -  X  ,  including  a  power  which  is  inadr  ss't.-le 

p 

since  the  integral  of  dfo/dt  (and  ’pas  well)  must  converge.  (In  fitting 
their  expressions  by  polynomials  of  the  form  of  Eq.  (3)  one  should  use 
only  the  range  of  energy  over  which  experimental  data  are  available. 

This  range  may  be  read  from  their  graphs.) 

Once  the  task  of  fitting  for  a  given  inner  shell  to  obtain  the 
constants  A.  has  been  accomplished,  one  can  obtain  the  predicted  cross 
section  for  ionization  of  that  shell  by  fast  electrons  by  carrying  out 
the  double  integral^ 

<(E,)-  J  ^  I 


dLt  cLr^ 


(7) 


where  0  =  E-|/E  and  E-j  is  the  energy  of  the  incident  electron.  Since  the 
doubly-differential  cross  section  d2a/dedn  is  related  to  the  GOS  through 

dzcr  _  2Tre4  i  df  }  (&) 

dLedV  E0£  £*1  ds 

Eqs.  (2)  and  (7)  lead  to 
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with  g(e)  given  by  Eq.  (3)  with  the  constants  A.  determined  from  experi¬ 
mental  data. 

The  integration  over  e  in  Eq.  (9)  may  be  carried  out  analytically. 
The  basic  integral  to  be  evaluated  may  be  written 


I„ui  -  J 


dx  x" 


(10) 
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Where  k  =  4  in  Eq.  (9),  but  could  be  taken  as  a  fitting  parameter  i- 
some  calculations.  Equation  (10)  can  also  be  expressed  as 


t  f  <4/  (a + 

(yz  +  sf 

where  A  = K ^z,  and 

T  r  d*  >* 

For  k  =  0,  1 ,  Eq.  (12)  gives 
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where  .  The  recursion  relation 

fh-sV 2  te-t 

T  — L )  r  T  4.  — J - i-i—  _  > 

3  1^-2  J 

may  be  used  to  find  any  of  the  J,,  once  dQ  arid  nave  been  found 
In  terms  of  the  integrals  I„  in  Eq.  (10),  Eq.  (9)  becomes 
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The  remaining  integration  over  momentum  transfer  ■  may  be  carried  out 
numerical ly. 

The  procedure  outlined  here  gives  a  reasonable  approach  for  obtaining 

GOS's  and  hence  inner-shell  ionization  cross  sections,  from  the  relatively 

2  3 

abundant  photoionization  data  ’  covering  a  wide  range  of  elements.  There 
remains,  however,  the  task  of  comparing  the  results  obtained  by  this 
procedure  with  differential  inverse  mean  free  paths  calculated  using 
theoretical  GOS  values,  e.g.,  for  aluminum^,  and  experimental  data  on  inner- 

5 

shell  ionization  by  electrons. 
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